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DESIGN OF POLYKETIDE SYNTHASE GENES 

FIELD OF INVENTION 

The present invention provides methods for the analysis of polyketides and the design of 
polyketide synthase genes. The invention relates to the fields of computational analysis, 
chemistry, molecular biology, and medicine. 

BACKGROUND OF THE INVENTION 

The class of compounds known as polyketides is a large family of diverse compounds 
synthesized primarily from 2-carbon unit building block compounds through a series of 
condensations and subsequent modifications. Polyketides occur in many types of organisms, 
including fungi and mycelial bacteria such as the actinomycetes. There are a wide variety of 
polyketide structures, and the class of polyketides encompasses numerous compounds with 
diverse activities. Epothilone, erythromycin, FK-506, FK-520, megalomicin, narbomycin, 
oleandomycin, picromycin, rapamycin, spinocyn, and tylosin are examples of such compounds. 

Given the difficulty in producing polyketide compounds by traditional chemical 
methodology, and the typically low production of polyketides in wild type cells, there as been 
considerable interest in finding improved or alternate means to produce polyketide compounds. 
See PCT Publication Nos. WO 95/08548; WO 96/40968; WO 97/02358; and 98/27203; Unites 
States Patent Nos. 5,962,290; 5,672,491; and 5,712,146; Fu et aL, Biochemistry 33: 9321-9326 
(1994); McDaniel etal, Science 262:1546-1555 (1993); and Rohr, Angew. Chem. Int. Ed. Engl. 
34(8): 881-888 (1995), each of which is incorporated herein by reference. 

Polyketides are synthesized in nature by polyketide synthase (PKS) enzymes. These 
enzymes, which are complexes of multiple large proteins, are similar to the synthases that 
catalyze condensation of 2-carbon unit building block compounds in the biosynthesis of fatty 
acids. The genes that encode PKS enzymes usually consist of three or more open reading 
frames (ORFs). Two major types of PKS enzymes are known that differ in their composition 
and mode of synthesis. These two major types of PKS enzymes are commonly referred to as 
Type I or "modular" and Type n or "iterative" PKS enzymes. 
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Modular PKSs produce many different polyketides, including a large number of 12-, 14-, 
and 16-membered macrolide antibiotics including erythromycin, megalomicin, methymycin, 
narbomycin, oleandomycin, picromycin, and tylosin. Each ORF of a modular PKS can 
comprise one, two, or more "modules" of ketosynthase activity, each module of which consists 
of at least two (if a loading module) and more typically three (for the simplest extender module) 
or more enzymatic activities or "domains." These large multifunctional enzymes (>300,000 
kDa) catalyze the biosynthesis of polyketide macrolactones through multistep pathways 
involving decarboxylative condensations between acyl thioesters followed by cycles of varying 
/J-carbon processing activities (see O'Hagan, D., The polyketide metabolites, E. Horwood, New 
York, 1991, which is incorporated herein by reference). 

During the past half decade, the study of modular PKS function and specificity has been 
greatly facilitated by the plasmid-based Streptomyces coelicolor expression system developed 
with the 6-deoxyerythronoiide B (6-dEB) synthase (DEBS) genes (see Kao et al., Science, 265: 
509-512 (1994), McDaniel etal 9 Science 262: 1546-1557 (1993), and U.S. Patent Nos. 
5,672,491 and 5,712,146, each of which is incorporated herein by reference). The advantages to 
this plasmid-based genetic system for DEBS are that it overcomes the tedious and limited 
techniques for manipulating the natural DEBS host organism, Saccharopolyspora erythraea, 
allows more facile construction of recombinant PKSs, and reduces the complexity of PKS 
analysis by providing a "clean" host background. This system also expedited construction of a 
combinatorial modular polyketide library in Streptomyces (see PCT publication No. WO 
98/49315, incorporated herein by reference). 

The ability to control aspects of polyketide biosynthesis, such as monomer selection and 
degree of 0-carbon processing, by genetic manipulation of PKSs has stimulated great interest in 
the combinatorial engineering of novel antibiotics (see Hutchinson, Curr. Opin. Microbiol 1: 
319-329 (1998); Carreras and Santi, Curr. Opin. Biotech. 9: 403-41 1 (1998); and U.S. Patent 
Nos. 5,962,290; 5,712,146; and 5,672,491, each of which is incorporated herein by reference). 
This interest has resulted in the cloning, analysis, and manipulation by recombinant DNA 
technology of genes that encode PKS enzymes. The resulting technology allows one to 
manipulate a known PKS gene cluster either to produce the polyketide synthesized by that PKS 
at higher levels than occur in nature or in hosts that otherwise do not produce the polyketide. 
The technology also allows one to produce molecules that are structurally related to, but distinct 
from, the polyketides produced from known PKS gene clusters. 
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Polyketides are assembled by polyketide synthases through successive condensations of 
activated coenzyme-A thioester monomers derived from small organic acids such as acetate, 
propionate, and butyrate. Active sites required for condensation include an acyltransferase 
(AT), acyl carrier protein (ACP), and beta-ketoacylsynthase (KS). Each condensation cycle 
results in a |3-keto group that undergoes all, some, or none of a series of processing activities. 
Active sites that perform these reactions include a ketoreductase (KR), dehydratase (DH), and 
enoylreductase (ER). Thus, the absence of any beta-keto processing domain results in the 
presence of a ketone, a KR alone gives rise to a hydroxyl, a KR and DH result in an alkene, 
while a KR, DH, and ER combination leads to complete reduction to an alkane. After assembly 
of the polyketide chain, the molecule typically undergoes cyclization(s) and post-PKS 
modification (e.g. glycosylation, oxidation, acylation) to achieve the final active compound. 

To illustrate the synthesis of a macrolide by a modular PKS (see Cane et aL, Science 
282: 63 (1998), incorporated herein by reference), one can refer to the PKS that produces the 
erythromycin polyketide (6-deoxyerythronolide B synthase or DEBS; see U.S. Patent No. 
5,824,513, incorporated herein by reference). In the modular DEBS PKS enzyme, the enzymatic 
steps for each round of condensation and reduction are encoded within a single "module" of the 
polypeptide (i.e., one distinct module for every condensation cycle). As shown in Figure 1, 
DEBS consists of a loading module and 6 extender modules and a chain terminating thioesterase 
(TE) domain within three extremely large polypeptides encoded by three open reading frames 
(ORFs, designated eryAI, eryAII, and eryAIII). 

Each of the three polypeptide subunits of DEBS (DEB1 , DEBS2, and DEBS3 in 
Figure 1) contains 2 extender modules. DEBS1 additionally contains the loading module, and 
DEBS3 contains the TE domain. Collectively, these proteins catalyze the condensation and 
appropriate reduction of one propionyl CoA starter unit and six methylmalonyl CoA extender 
units. Modules 1, 2, 5, and 6 contain KR domains; module 4 contains a complete set, 
KR/DH/ER, of reductive and dehydratase domains; and module 3 contains no functional 
reductive domain. Following the condensation and appropriate dehydration and reduction 
reactions, the enzyme bound intermediate is lactonized by the TE at the end of extender module 
6 to form 6-dEB (compound 1 in Figure 1). 



More particularly, the loading module of DEBS consists of two domains, an acyl- 
transferase (AT) domain and an acyl carrier protein (ACP) domain. In other PKS enzymes, the 
loading module is not composed of an AT and an ACP but instead utilizes a partially inactivated 
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KS, an AT, and an ACP. This partially inactivated KS is in most instances called KS Q , where 
the superscript letter is the abbreviation for the amino acid, glutamine, that is present instead of a 
cysteine in the active site that is believed to be required for condensation activity. Although the 
KS Q domain lacks condensation activity, it retains decarboxylase activity. The AT domain of 
the loading module recognizes a particular acyl-CoA (propionyl for DEBS, which can also 
accept acetyl) and transfers it as a thiol ester to the ACP of the loading module. Concurrently, 
the AT on each of the extender modules recognizes a particular extender-CoA (methylmalonyl 
for DEBS) and transfers it to the ACP of that module to form a thioester. Once the PKS is 
primed with acyl- and malonyl-ACPs, the acyl group of the loading module migrates to form a 
thiol ester (trans-esterification) at the KS of the first extender module; at this stage, extender 
module 1 possesses an acyl-KS and a methylmalonyl ACP. The acyl group derived from the 
loading module is then covalently attached to the alpha-carbon of the malonyl group to form a 
carbon-carbon bond, driven by concomitant decarboxylation, and generating a new acyl- ACP 
that has a backbone two carbons longer than the loading unit (elongation or extension). The 
growing polyketide chain (various intermediates are shown in Figure 1) is transferred from the 
ACP to the KS of the next module, and the process continues. 

The polyketide chain, growing by two carbons each module, is sequentially passed as a 
covalently bound thiol ester from module to module, in an assembly line-like process. The 
carbon chain produced by this process alone would possess a ketone at every other carbon atom, 
producing a polyketone, from which the name polyketide arises. Commonly, however, 
additional enzymatic activities modify the beta keto group of the polyketide chain to which the 
two carbon unit has been added before the chain is transferred to the next module. Modules may 
contain additional enzymatic activities as well, such as methyl transferase domains, but there are 
no such additional activities in DEBS. 

Once a polyketide chain traverses the final extender module of a modular PKS, it 
encounters the releasing domain or thioesterase found at the carboxyl end of most PKSs. Here, 
the polyketide is cleaved from the enzyme and cyclyzed. The resulting polyketide can be 
modified further by tailoring or modification enzymes; these enzymes add carbohydrate groups 
or methyl groups, or make other modifications, i.e., oxidation or reduction, on the polyketide 
core molecule. For example, the final steps in conversion of 6-dEB to erythromycin A include 
the actions of a number of modification enzymes, such as:' C-6 hydroxylation, attachment of 
mycarose and desosamine sugars, C-12 hydroxylation (which produces erythromycin C), and 
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conversion of mycarose to cladinose via O-methylation. These modifications in various 
combinations result in erythromycins A (compound 2 in Figure 1), B, C, and D. 

While the detailed understanding of the mechanisms by which PKS enzymes function 
and the development of methods for manipulating PKS genes have facilitated the creation of 
novel polyketides, there remain substantial impediments to the creation of novel polyketides by 
genetic engineering. One such impediment is the availability of PKS genes. Many polyketides 
are known but only a relatively small portion of the corresponding PKS genes have been cloned 
and are available for manipulation. Moreover, in many instances the producing organism for an 
interesting polyketide is obtainable only with great difficulty and expense, and techniques for its 
growth in the laboratory and production of the polyketide it produces are unknown or difficult or 
time-consuming to practice. Also, even if the PKS genes for a desired polyketide have been 
cloned, those genes may not serve to drive the level of production desired in a particular host 
cell. 

If there were a method to produce a desired polyketide without having to access the 
genes that encode the PKS that produces the polyketide, then many of these difficulties could be 
ameliorated or avoided altogether. The present invention meets this need. 

SUMMARY OF THE INVENTION 

In one embodiment, the present invention provides methods for the computational 
analysis of polyketides and the computer-assisted design of PKS genes. 

In a first aspect, the present invention provides a method for representing the structure of 
a polyketide and/or a PKS gene that encodes the PKS that produces the polyketide by 
alphanumeric symbols that facilitates computer assisted analysis. 

In a second aspect, the present invention provides a database of polyketides and 
corresponding PKS genes that can be rapidly searched and information extracted for a variety of 
applications. More particularly, this database can include, in one mode, all known polyketides; 
and in another mode, the polyketides, optionally including all intermediates, produced by all 
known PKS genes or a subset thereof. 
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In a third aspect, the present invention provides a method for predicting the structure of a 
PKS and its corresponding genes from the structure of a polyketide. 

In a fourth aspect, the present invention provides a method for designing novel PKS 
genes capable of producing a desired polyketide. This aspect of the invention is directed to the 
design and specification of PKS genes via the recombining of modules or portions of modules or 
sets of modules from already known and available PKS genes. In one mode, all possible PKS 
genes encoding a desired polyketide from a set of genes in a database are generated. In another 
mode, only a subset of such possible PKS genes is generated based on one or more parameters 
selected by the user. More particularly, a rating system is provided to sort the PKS genes 
designed for a particular target polyketide based on any one or more of several criteria, including 
number of non-native module interfaces, number of non-native protein interfaces, and other 
parameters as more particularly described below or selected by the user. 

In another embodiment, the present invention provides methods and reagents for 
preparing novel PKS genes that encode PKS enzymes that produce a desired polyketide. 

In a first aspect, the present invention provides a library of recombinant DNA 
compounds, wherein each member of said library encodes a module of a PKS or portions of 
modules or sets of modules having a desired specificity, and the library as a whole encompasses 
all of the members of a desired class of specificities. 

In a second aspect, the present invention provides a method for assembling a PKS gene 
cluster that encodes a PKS that produces a desired polyketide from known and available PKS 
genes other than the naturally occurring PKS genes that produce the polyketide in nature. 

These and other embodiments, modes, and aspects of the invention are described in more 
detail in the following description, the examples, and claims set forth below. 

BRIEF DESCRIPTION OF THE FIGURES 

Figure 1 shows a schematic representation of the PKS enzyme that synthesizes 6- 
deoxyerythronolide B (6-dEB, compound 1). The PKS is composed of three proteins, DEBS1, 
DEBS2, and DEBS3, each of which is represented by an arrow and contains two or more 
modules. Each module is represented by a solid line, and the domains in each module are shown 
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inside the arrow. Various intermediates produced during the synthesis are also shown, as are the 
structures of erythromycins A (compound 2), B, and D resulting from modification of 6-dEB. 

Figure 2 shows an illustrative set of 2-carbon unit monomers present in macrocyclic 
polyketides; these monomers can be used to represent polyketide backbone diversity generated 
by commonly used starter and extender units (malonyl CoA and methylmalonyl Co A) and the 
condensation and reduction reactions mediated by PKS enzymes. 

Figure 3 shows a representation of 6~dEB by molecular graph, CHUCKLES notation, 
and SMILES notation. The CHUCKLES notation uses the 2-carbon unit monomers shown in 
Figure 2. In the CHUCKLES notation, the order of attachment of monomers is designated by 
the order in which monomers are listed, and the attachment points within the monomers are 
specified in their definitions. In the SMILES notation, adjacent monomers are attached via 
single (covalent) bonds depicted by dashes. The cyclization bond is represented by the index 1 
adjacent to the Start and Close monomers. 

Figure 4 is a flowchart and block flow diagram in five parts designated A-E, inclusive. 

Flowchart Figure 4A is a block flow diagram of a computer system to design a novel 
PKS (and corresponding genes). 

Flowchart Figure 4B is a block flow diagram wherein the "Computer Program'* block (2) 
of Flowchart Figure 4A is further defined. 

Flowchart Figure 4C is a block flow diagram wherein the "Design novel hybrid PKS 
genes from library for TARGET" block of Flowchart Figure 4B is further defined. 

Flowchart Figure 4D is a block flow diagram wherein the "align TARGET with 
STARTER; copy to ALIGNMENT" block of Flowchart Figure 4C is further defined. 

Flowchart Figure 4E is a block flow diagram wherein the "Rate novel hybrid designs" 
block (3) of Flowchart Figure 4B is further defined. 

Figure 5 shows a flowchart of a matching method for the generation of the CHUCKLES 
strings used for all polyketides in a library. 



WO 01/92991 



8 



PCT/US01/17352 



DETAILED DESCRIPTION OF THE INVENTION 



Because polyketides synthesized by modular PKS genes are built by the enzymatically 
controlled addition of primarily 2-carbon unit monomers and, to a lesser extent, other more 
complex monomers, each polyketide may be represented as a string of 2-carbon unit and other 
monomers. These monomers represent the portion of the polyketide backbone structure as a 
result of the incorporation of various starter and extender units (malonate, methyl malonate, etc.) 
and the subsequent chemical reactions. 

These reactions include: 

(1) condensation reactions, of which there are three basic reactions: malonyl-CoA 
condensation arid methylmalonyl-CoA condensation with the branched methyl having either R 
or S stereochemistry; and 

(2) reduction reactions, of which there are five basic reactions: no reduction (ketone 
preserved), keto-reduction (to yield a hydroxyl having either R or S stereochemistry), 
dehydration (trans double bond), and enoyl-reduction (to yield a methylene). 

An illustrative set of the basic monomers that can be used to represent a polyketide 
structure (and their corresponding symbols) comprises: 

OH OH OH OH 





O 





and 
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A miscellaneous monomer, Q, can be used to denote a portion of the polyketide structure that 
cannot be assigned by monomers A-N. 

The monomer set shown above and in Figure 2 does not represent the actual monomers 
incorporated during biosynthesis. Instead, these monomers include a carbon from two different 
biosynthetic monomers. This is best explained using a polyketide fragment depicted below. 




CH3 CH3 

The fragment includes two two-carbon units, i and i+1 and part of a third two-carbon unit, i-1 
that were incorporated into the polyketide during biosynthesis. The i-th extender module 
attaches the two carbon biosynthetic unit whose backbone carbons are designated as alphas and 
beta* and the second extender module attaches the two carbon biosynthetic unit whose backbone 
carbons are designated as alphaj+i and beta,+i. Using the monomer set shown above, this 
fragment consists of monomer A (derived from the beta carbon added in module i+1 and the 
alpha carbon added in module i) and another monomer A (derived from the beta carbon added in 
module i and the alpha carbon added in module i+1). 




CH3 CH3 

1 I 1 I 

A A 

The fifth carbon designated beta H1 remains unassigned and will depend on the identity of the 
two-carbon biosynthetic unit that is incorporated in the polyketide by module i+2. 

The set of monomers shown in Figure 2 can be expanded to include other starter and 
extender units, of which there are many. Such starter and extender units include, for example 
but without limitation, hydroxymalonate (e.g., niddamycin), methoxymalonate (e.g. FK-520), 
ethylmalonate (e.g., FK-520), amino acids or amino acid derivatives that are incorporated into 
polyketides by the action of a non-ribosomal peptide synthase (e.g., thiazole in epothilone and 
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pipecolate in rapamycin), or other units incorporated by, for example, an AMP ligase (e.g., the 
dihydroxycylohexyl moiety in rapamycin, FK-506, and FK-520) or a soluble CoA ligase. An 
illustrative set of additional starter and extender units includes: 




=J ' ; = K ';ami /^y' = M ' 

R A R 

where R can be anything other than hydrogen or methyl (e.g., allyl, butyl, ethyl, hexyl, hydroxyl, 
isobutyl, and methoxy). 

The set of monomers can also include post-PKS modifications, such as hydroxylation, 
methylation, epoxidation, glycosylation, or addition of intra-macrocyclic fused rings making the 
system polycyclic. Also, a variety of methods are known for the incorporation of unusual starter 
and or extender units in polyketide synthases (see, e.g., PCT Publication Nos. WO 97/02358; 
WO 99/03986; WO 98/01546; and WO 98/01571, each of which is incorporated herein by 
reference, and the monomer set can include such units. 

By viewing polyketides as composed of sets of distinct monomers, one can in 
accordance with the present invention define a polyketide as a string of alpha-numeric symbols 
to facilitate computer analysis. In one method, a modified CHUCKLES methodology for 
representing polyketides is used. The CHUCKLES methodology (see Siani et ai, 
"CHUCKLES: a method for representing and searching peptide and peptoid sequence," J. 
Chem. Inf. Sci. 34: 588-593 (1994) which is incorporated herein by reference) for representing 
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peptides and related oligomers allows monomers to be strung together such that the molecular 
graph for the basic macrocycle can be generated from the string of monomers. 

For example, using the set of monomers comprising A-N described above, the 
erythromycin macrocycle or 6-dEB can be represented as ADGJDD. This string of 
alphanumeric symbols is also referred to as the CHUCKLES string. Figure 3 depicts the 
relationship between the CHUCKLES string, the SMILES string, and the actual molecular 
structure of 6-dEB. The CHUCKLES string for 6-dEB can be annotated to represent the 
structure of erythromycin A: A(l-lactone closure,2-hydroxyl)DGJ(2-hydroxyl)D(l-glycosyl) D(I- 
glycosyl). Thus, ring closure (cyclization) and post-synthetic modifications (glycosylation and 
hydroxylation), and non-standard units where applicable (there are none in 6-dEB and 
erythromycin) are entered between parentheses after each monomer. Another example is an 
annotated CHUCKLES string for epothilone B: ME(l-lactone-closure)M(epoxide)LJDG(2- 
methylation)E. As above, cyclization, post-synthetic modifications (epoxide formation), and 
non-standard units (methyl at C-4) are entered between parentheses after each monomer. 

In another aspect of the present invention, a database of polyketides is provided. In one 
aspect of the present invention, the polyketides are represented by a string of defined monomers. 
In one embodiment, the monomers are selected from a group consisting of: 
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In another embodiment, polyketides are represented by the monomers A-N as well as 
additional monomers selected from the group consisting of 




= J ' ; / Y = K ' ; = M ' 

R R & 

where R can be anything other than hydrogen or methyl. 

The string of monomers can be represented as a linearized structure or as a string of 
symbols. For example, the erythromycin can be represented as its aglycone, 6-dEB, as 




or as a string of symbols, ADGJDD. Optionally, the string of symbols can be annotated as "A(l- 
lactone closure,2-hydroxyl)DGJ(2-hydroxyl)D(l-glycosyl) D(I-glycosyl)" to more fully capture 
the erythromycin structure. This set of annotated strings is referred to as a "coded library" or a 
"coded" database of the present invention. 

In an illustrative embodiment, the polyketide database consists of the polyketides 
described in current literature (Journal of Antibiotics (1981-present), Journal of Natural 
Products) and various databases (Chemical Abstracts CAPlus, AntiBase). All unique 
macrocyclic polyketides are converted to the modified CHUCKLES format. Of the -1000 novel 
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polyketides obtained, only -200 different strings of monomers and unique macrocycles are 
needed to represent the much larger collection of polyketides in the database, because many of 
the differences between the naturally-occurring polyketides are due to different glycosyl (sugar) 
groups attached at different positions on the macrocycle. 

Thus, a macrocyclic polyketide can be converted to a string of 2-carbon monomers by 
mapping the monomers onto the polyketide. This can be performed manually or with computer 
assistance. First, any sugar moieties are conceptually removed by hydrolysis and any lactones 
(bond between the ketone and oxygen) are hydrolyzed thus generating a linearized structure of 
the backbone of the polyketide. Generally, this leaves a carboxy carbon at one end of the linear 
molecule and a hydroxyl at the other. The polyketide is then "sequenced" manually or in silico 
from the end containing the carboxy carbon, the end corresponding to the last monomer added 
by the PKS before synthesis is complete. This end serves as a convenient handle from which to 
start the mapping process. Although closing of the lactone often occurs between the two ends of 
the polyketide, this is not always the case. However, the last ketone added by the PKS is almost 
always involved in macrolactone formation and so serves as a more convenient handle than the 
hydroxyl for commencing sequencing. 

The manual or in silico sequencing is performed by matching the monomers, one at a 
time, while traversing the macrocyclic backbone. First the carboxy carbon is skipped, and an 
attempt is made to match each of the monomers in the monomer set selected (i.e., monomer set 
A-N in Figure 2) against the next two carbons in the macrocycle. The match takes into account 
carbon, oxygen, and no substitution at each backbone position, chirality at each backbone 
position, and bond order between the two backbone carbons. 

If the sequencing is performed in silico, the method is referred to as back-translation and 
involves converting a molecular graph into a string of monomers. First, the monomer library is 
converted to SMARTS format. SMARTS is a superset of the SMILES language that specifies a 
pattern in a molecular graph (Daylight Software Manual: Theory; Daylight Chemical 
Information Systems; Irvine, CA 1993, incorporated herein by reference). SMARTS permits 
one to specify a variable number or a limit on the number of covalent bonds to non-hydrogen 
atoms from a particular atom. In contrast, SMILES assumes that the unspecified valences are 
hydrogens. For example, the SMILES string for monomer A is [C@@H](0)[C@H](C). The 
oxygen may be bonded to any other single atom; if the atom is not specified, it is assumed be a 
hydrogen. In the SMARTS string for monomer A, [C@@H](0;D2])[C@H]([CH3]), one can 
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specify the exact number of hydrogens on some atoms (e.g., "CH3"). In addition, the "[0;D2]" 
indicates the oxygen is bonded to two (from D2) non-hydrogen atoms, in this case the first 
carbon and some other unspecified atom. This allows matching and distinction of post- 
modification moieties attached to the oxygen as well as additional cyclizations (six member 
rings can occur within the macrocycle; e.g., rapamcyin). Thus, the SMARTS notation allows 
pattern matching against the polyketide molecular graph. 

When a match occurs, the atoms that match are tagged as part of a superset and labeled 
with the monomer name. Any atoms that are connected to the monomer that are not part of the 
macrocycle are tagged for identification as special precursor units (e.g., ethylmalonate instead of 
methyl malonate or malonate), or post-synthetic modification moieties (e.g., sugars, CCHO, 
hydroxylation, methylation). If all the atoms and boncls of the monomer cannot be identified, 
the monomer is given a designation to indicate the lack of identification (e.g., Q for question 
mark). These Q monomers can be used to identify monomers that are the site of post-PKS 
modifications that mask the function of the PKS module that generated that portion of the 
polyketide or that are not in the monomer set and so prevent the correlation of a particular 
segment of the backbone with one of the monomers in the monomer set. 

After a particular 2-carbon unit is identified, the next two carbons are processed the same 
way. This is repeated until all the backbone carbons are identified and labeled as monomers. 
When all two-carbon units are identified, one has generated an ordered sequence, or string, of 
monomers, which is a modified CHUCKLES string of the invention. Moieties corresponding to 
post-PKS modifications are appended to the monomer in the string as an annotation in 
parentheses. This method of sequencing maybe extended to include any type of monomer. 
Figure 5 shows a flow chart of this matching method for the generation of the CHUCKLES 
strings used for all polyketides in a library. 

The CHUCKLES string can be in the order corresponding to the direction of 
biosynthesis on the PKS or its reverse. Each CHUCKLES string has a one-to-many relationship 
with the PKS gene in the producing organism. Thus, while many different organisms can 
produce the same polyketide using the same or different PKS genes, each PKS gene generally 
produces only one PKS that produces only one polyketide (some AT domains can bind different 
CoAs, leading to the production of multiple polyketides from a single PKS). This allows one to 
design, from the polyketide structure, a set of PKS genes that would produce that polyketide. 
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Thus, the present invention provides methods and computational analysis tools for 
designing PKS genes to produce a desired polyketide. As an illustrative example, the present 
invention provides a computer program termed MORPH (see the Examples below) that can read 
the coded library (see the Examples below). An illustrative coded library consists of -200 
unique polyketide CHUCKLES strings. The user specifies the target polyketide, which is 
converted from molecular structure to a CHUCKLES string. 

The program then performs the following, starting with each library compound or string: 

(1) aligns library compound and target compound, emphasizing alignment of adjacent 
monomers common between the two; 

(2) fills in the gaps using all possible combinations from all library members; 

(3) counts number of non-natural inter-modular boundaries, 

(4) outputs all these alignments. 

The alignments are then sorted based on the number of non-natural inter-modular 
boundaries. 

This illustrative program allows one to design and find PKS genes that encode PKS 
enzymes that are combinations of two or more different PKS enzymes with the fewest inter- 
modular boundaries, and optionally the fewest inter-protein boundaries. Many other alternative 
embodiments are provided by the present invention. 

For example, one can include the naturally occurring PKS that produces the target 
polyketide in the coded library to allow components of that PKS to be incorporated into the 
design of a new PKS. Also, one can include in the coded library non-naturally occurring PKS 
enzymes, such as those produced and published in the scientific and patent literature to make 
novel polyketides, in the coded library. See, e.g., PCT publication Nos. WO 98/49315 and WO 
96/40968, both of which are incorporated herein by reference. 

This CHUCKLES-coded polyketide library can be stored in a computer file as a set of 
records. In one embodiment, each record contains the chemical name of the polyketide, the 
unannotated CHUCKLES (containing basic macrocyclic monomers), the annotated 
CHUCKLES (containing basic macrocyclic monomers with information about post-PKS 
modifications), the producing organism(s), and other information (e.g., linearized representation 
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of the polyketide structure, the accession number of organisms or plasmids that have been 
deposited, gene sequence information, and references). 

The MORPH program can read in the polyketide library entries to an array or list of dat 
structures, where each entry data structure contains all or a selected subset of the fields in each 
library record. The MORPH program then reads in the CHUCKLES-coded TARGET 
polyketide from the user. This TARGET may optionally be blocked from the library so that it : 
not used as a STARTER or left in the library, i.e., if it is only distantly related to other known 
polyketides, or some modules could be useful in designing novel PKS genes, or it is desired to 
replace only certain PKS modules. This program could also be used for analoging at a particul 
position via wild-cards defined as part of the TARGET sequence by the user. 

Each member of the coded PKS library can be selected as a STARTER unit. Thus, 
during a rim, all library members can be given an equal chance as STARTER units. After a 
STARTER is chosen, the TARGET is aligned with it. See Flowchart Figure 4D. Any method 
of alignment can be used such that the maximal number of adjacent STARTER modules is use< 
in the final alignment. After the maximal adjacent modules are used in the ALIGNMENT, 
smaller adjacent sets or individual modules from the STARTER are used to fill in the gaps: 
There may be several alignments that are equally good based on the attempt to optimize the 
number of adjacent modules. For example, if the TARGET contains the "JDG" substring, then 
6-dEB, identified with the A1D2G3 J4D5D6 CHUCKLES string, may align as 



TARGET 


JDG 


6-dEB 


J4D2G3 


TARGET 


JDG 


6-dEB 


J4D5G3. 



Both of these alignments have different maximal adjacent modules, with the 
same length of two (D2G3 in the first and J4D5 in the second). Accordingly, either 
alignment could be used as STARTERs. 

With the optimized alignment from the STARTER, other library entries are 
systematically used to complete the alignment, or fill in the gaps. This part may be performed 
on either the optimized ALIGNMENT described above, or the ALIGNMENT without the sing] 
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modules from the STARTER; the removal of the individual modules opens up more space into 
which larger pieces of the FILLER might be placed. The first library entry is designated as the 
FILLER. If the FILLER is the same as the STARTER, the next library entry is used as the 
FILLER. This library entry is flagged as the CURRENT_FILLER_LIBRARY_ENTRY. The 
same method for finding maximally adjacent modules and then smaller sets or single modules is 
used to fill the gaps in ALIGNMENT from the FILLER. If not all the gaps are filled in the 
ALIGNMENT, then the next library entry is used as a new source; that is, it is designated as the 
FILLER, and the gaps are filled further. This is repeated until the ALIGNMENT is complete or 
the end of the library is reached. 

Assuming all modules in the TARGET are represented in the library, the ALIGNMENT 
is eventually completely filled. The completed alignment is then written to an output file on the 
computer disk. When the ALIGNMENT is complete, or there are no more FILLERS in the 
library, the TARGET and STARTER alignment are re-copied to ALIGNMENT. The 
CURRENT_FILLER_LIBRARY_ENTRY is incremented, and a new attempt to fill in the gaps 
is started. 

When the CURRENT_FILLER_LIBRARY_ENTRY has reached the end of the library, 
the ALIGNMENT is wiped, and a new STARTER is chosen. The above process is then 
repeated for the next STARTER. When all library entries have been used as starters, then all 
feasible novel polyketide synthases have been generated and written to the computer file. The 
novel PKSs are then read back into memory and can be further evaluated. An illustrative 
evaluation process involves: 

(1) counting the non-native inter-module interfaces, and 

(2) counting the number of native inter-protein interfaces (for known and 
annotated gene sequences). 

The novel PKSs are then sorted based on these two numbers, giving higher priority to the non- 
native inter-module interfaces. In this mode, the goal is to identify those novel PKSs that contain 
the fewest non-native interfaces. 

By providing methods and means for the computer-assisted analysis of polyketides and 
PKS genes, the present invention greatly facilitates the identification and production of new 
polyketides with useful activities. Those of skill in the art will appreciate that while the 
invention is in part illustrated in the Examples below with respect to the design of new PKS 
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genes for known polyketides, the invention can also be used to design PKS genes for novel 
polyketides. In this embodiment, one simply provides the structure of the novel polyketide to 
the MORPH or other program of the invention to generate the desired PKS genes. 

Moreover, while the invention is exemplified below by designing new PKS genes 
composed of the coding sequences for one or more complete modules of two or more different 
PKS genes, partial modules can also be employed. With the appropriate choice of monomer sets 
and corresponding coding of the library to be searched, one can generate new PKS gene designs 
that take advantage of the potential to fuse one PKS gene coding sequence to another at a site 
corresponding to an intra-modular junction. In another embodiment, one can use "wild-cards" 
in the encoded polyketide or library to take advantage of known or predicted SAR. Thus, if one 
knows that a particular position in a polyketide can be varied (i.e., a hydrogen, methyl, or ethyl 
group at a location determined by an AT domain of a particular module, or a hydroxyl or keto 
group at a location determined by the presence or absence of a KR domain in a particular 
module) then one can use a wild-card monomer designation in the polyketide CHUCKLES 
string to generate PKS genes that produce each of the desired variants. 

The methods of the invention have diverse application in addition to the design of new 
PKS genes. As but one illustrative example, the methods of the invention can be used to design 
methods to produce a desired compound. Organic molecules containing stereochemical centers 
are useful for a number of purposes, including use as synthetic or semi-synthetic intermediates. 
The preparation of such intermediates by organic synthesis can be extremely time consuming 
and expensive. An alternative source of such intermediates is via specific degradation of a 
polyketide, and the present invention provides computer-assisted means for designing such 
production methods. 

Thus, certain functional groups of polyketides are susceptible to bond cleavage by 
specific chemical reactions that do not affect other functional groups. For example, carbon- 
carbon double bonds can be specifically cleaved by permanganate without affecting other 
functional groups normally in polyketides, such as ketones, alcohols, and lactones. Likewise, 
the Baeyer- Villager reaction converts a ketone to an ester (lactone) without affecting other 
groups of the aglycone. In accordance with the methods of the invention, one can assemble a 
library of polyketides in a database that can be addressed with a query describing a particular 
chemical reaction to generate all of the degradation products produced by that reaction upon 
each of the polyketides in the library. The degradation fragments thus generated serve as a 
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library of the invention that can be sorted by properties, such as size, number and type of 
stereochemical centers, functional groups, or other factors, and searched for useful compounds. 
Moreover, the functional groups on the ends of the fragments generated (or at other locations) 
can also be converted to other functional groups by chemical reactions (optionally employing 
protecting groups on other functional groups), and the database of compounds can be expanded 
to include the compounds produced by such reactions. 

From even a modest library of -200 compounds, one can in this manner generate using 
the methods of the invention, two to three times as many valuable chemical intermediates. Once 
such an intermediate is identified, the organism that produces the polyketide from which the 
fragment is derived is fermented, the polyketide isolated in bulk, the chemical reaction 
performed, and the desired degradation produces) isolated and used. In this manner, the present 
invention makes available a wide variety of useful products otherwise unattainable. 

Thus, the present invention has wide application in the fields of chemistry, particularly 
medicinal chemistry, molecular biology, and medicine. Those of skill in the art will recognize 
these and other benefits and applications provided by the present invention. Thus, the following 
examples are given for the purpose of illustrating the present invention and shall not be 
construed as being a limitation on the scope of the invention or claims. 

EXAMPLE 1 

The MORPH Program 

This example provides the source code for an illustrative MORPH program of the 
invention. The MORPH program is a command line driven program that runs on a UNIX 
system. The program can be run from a shell script, such that the user fills in the entire 
command ahead of time, then post-processes the output file with UNIX utilities including sort, 
egrep, and uniq. 

The command line appears as follows: 

motph3 -1 libraryfile -n targetname -ttargetsequence [-x X-wildcards] [-y Y-wildcards] [-z Z- 
wildcards]. 
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The library file is the name of the text file described below in Example 2. The target 
name is a user-defined identifier to distinguish this target from the library members (e.g., 
epothiloneD). The target sequence is a string of monomers that represent the CHUCKLES- 
encoded target polyketide (e.g., MEMLJDGE). Generally, if the target sequence is in the 
library, it is commented out from the library so that the morph program does not find the target . 
itself. The three different wildcards, X, Y, and Z, are independent sets of monomers that can be 
included in the target sequence. 

The output from the morph program can be redirected to a file. This output file is then 
post-processed by (1) extracting the HIT lines with valid combinations of modules that yield the 
target, (2) sorting the.HITS based on alphanumeric content using the UNIX sort command, (3) 
running the UNIX uniq command which removes multiple copies of each HIT, leaving one copy 
of each, (4) sorting based on the number of pieces in the sequence of modules. Generally, the 
fewest number of pieces, which correspond to the fewest number of inter-modular interfaces, are 
desired; these will appear at the top of the output. 

Below are some illustrative examples of calls to the MORPH program from a shell script 
using epothilone as a target. The first example generates combinations that yield epothilone D: 

%morph3 -1 PKS.lib -n epoD -t MEMLJDGE > omorph3_epoD 

%egrep HIT omorph3_epoD | sort | uniq | sort +10 -1 1 > omorph3 epoD.uniq.sort 

The second example generates combinations that yield a derivative of epothilone D having a 

hydroxylatC-13: 

%morph3 -1 PKS.lib -n epoD-130H -t MEXLJDGE -x ABCD > oepoD-130H 
%egrep HIT oepoD-130H | sort | uniq | sort +10 -1 1 > oepoD-130H.uniq.sort 

The third example generates an epothilone having wildcards (set 1): 

%morph3 -1 PKS.lib -n epoD-setl -t MEXYZDgE -x ABCD -y LEFIN -z JACGM > 

oepoD-setl 

%grep HIT oepoD-setl | sort | uniq ) sort +10 -1 1 > oepoD-setl.uniq.sort 
The fourth example generates an epothilone having another set of wildcards (set 2): 

%morph3 -1 PKS.lib -n epoD-set2 -t MEXYZDgE -x JK -y EF -z JACGM > oepoD-set2 
%grep HIT oepoD-set2 | sort ) uniq | sort +10 -11 > oepoD-set2.uniq.sort 

MORPH in its current implementation operates at the monomer level and thus does not 
handle intra-modular modifications/splitting. Future implementations could convert the 
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CHUCKLES-encoded strings into the corresponding and equivalent SMILES and then perform 
more complex chemical analysis of the PKS molecular graphs. Currently, inter-modular double 
bonds are present in the library, but are ignored by the program. These bonds can introduced 
post-biosynthetically and the exact source is generally unknown. 

The source code for MORPH is found in Appendix A (version 3.0) and B (version 4.0) 
(deposited in the microfiche appendix). 



EXAMPLE 2 

Illustrative Polyketide Library 

This example provides the contents of an illustrative CHUCLKES encoded polyketide 
library. The first column provides the name of the polyketide; the second the CHUCKLES 
string; the third the annotated CHUCKLES string; and the fourth the source organism. Entries 
under annotated CHUCKLES and source organism are not complete for all of the polyketides. 
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POLYKETIDE 


CHUCKLES 


annotated-CHUCKLES 


SOURCE 
ORGANISM 


3-acety]-4"- 
butyltylosin 


FMNGODF 






#aculeximycin 


RRNRSMRSRSS 
SSRLSSN 


RR(2-ethyl)NRSMRS( 1 - 

glycosyl)RSSS(2-hydroxyl)SR(l- 

glycosyl)LSSN(2-ethyl) 




albocycline-Ml- 

ingramycin- 

TA2407- 

cineromycinB- 

U28010-SR2077 


BLME=JN 


tvw n t r i /-« \t-i/i j i -\ t/ci 

BLM( 1 >2-epoxy)E( 1 -methoxy)= J(2- 
hydroxyl)N 




albocycline-M2 


BLME=JN 


B(2-hydroxyl)LME(l -methoxy)=J(2- 
hydroxyl)N 




albocycline-M3 


BSME=JL 


BSME( 1 -methoxy)=J(2-hydroxy 1)L 




albocycKne-M5 


BSME=JN 


BSME(1 -methoxy)=J(2-hydroxyi)N 




albocycline-M6 


BLME=JL 


BL(2-hydroxyl)ME(l -methoxy)=J(2- 
hydroxyl)L 




albocycline-M7 


BLME=JN 


BL(2-hydroxyl)ME( 1 -methoxy)=J(2- 
hydroxyl)N 




albocycline-M8 


BLME=Q 


BLME(l-methoxy)=Q 




aldgamycin 


BMLGJDL 


B(l-cyc)MLG(2-hydroxyl)JDL(2- 
cyc) 




amphotericinA 


CDNNLNNNNF 
CEFEALEE 


CDNNLNNNNF( 1 -glycosyl)C ( 1 -0- 
cyc,2-carboxylicacid)EF( 1 - 
cyc)EALEE 




#amphotericinB 


CDNNNNNNNF 
QQQEELEE 


CDNNNNNNNF(1 -glycosyl)C( 1 -0- 

cyc,2-carboxylicacid)EF(l- 

cyc)EALEE 




angiolam 


NMFJNSJIQLGA 
m 






aplyronineA 


BFJCENFFMEK 
AFNN 


B(1-C(=0)C)F(1- 
C(=0)C(C)N(C)C)JCENFFME(1- 
methoxy)KAF(l- 
C(=0)C(N(C)C)COC)NN 


sea hare Aplysia 
kurodai 


apoptolidin 


EFLMNAMMM 


EF(methoxy-l,hydroxy-2)LMNAMMM 


aurachinB 


MLMLM 


MLMLM 




aurachinC 


MLMLM 


MLMLM 




A59770 


QQKQQLJFCDN 


QQK(2-ethyl)Q=QU(2- 

hydroxyl)F(2-0-glycosyl)CD(2- 

hydroxyl)N 


Amycolatopsis 
orientalis 


A82548A- 
cytovaricin 


QQQKQNLJFDD 
N 






A83543A 


FLFQQQ 


FLFQQQ 


Saccharopolyspora 
spinosa 


AB023a 


NNNNNRSSSLS 
R 


NNNNNRSSSLSR 




AH -75 8 


RSNMURMN 


RS(2-methoxy)NMURMN(2- 
methoxy) 




bafilomcinD 


BNHCENMKCM 
N 


BNHCE(1 -macrocyc,2- 
methoxy)NMKCMN(2- 
methoxvWflceto-macrocyc) 


S. sp. 


bafilomycinAl 


BNHCENMKCM 
N 


BNHCE(l-macrocyc,2- 
methoxy)NMKCMN(2- 
methoxy)Q(keto-macrocyc) 
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POLYKETIDE 


CHUCKLES 


annotated-CHUCKLES 


SOURCE 

IJKijrAJNIJSM 


borrelidin 


SNMRUUUS 


SNM(2-N)RUUUS 




calyculin 


QFBBMmNm 


QF(1 -methoxy)BBMmNm 


Discodermia calyx 


candicidin- 

candeptin-ascosin- 

levorin-etc 


CDNNNNNNNF 
CEELDEE 


CDN^M^MNNFf 1 -glycosyl)C( 1 -0- 
cyc,2-carboxylicacid)EE(l -cyc)E(2- 
hydroxyl)LIEE 


S. griseus, S. 
canescus, S. 
levoris, S. 
viridoflavus, Stv. 
grisoviridum 


candidin 


QDNNNNNNNF 
CEFELIEE 


RDNNNNNNNF(l-glycosyl)C(l-0- 
cyc,2-carboxylicacid)EF( 1 -cyc)E(2- 
hydroxyl)LBEE 




carbomycin 


ENNCODF 


EN( 1 ,2-epoxy)NHOF( 1 -glycosyl,2- 
methoxy)F(l-C(=0)C) 




carbomycinB- 
magnamycinB 


FNNGOFF 


FNNGOF(l-glycosyl,2-methoxy)F(l- 
C(=0)C) 




carbomycin-A- 

magnamycin- 

deltamycmA4- 

NSC51001-PS97- 

3628-WC3628 


ENNHOFF 


ENNHO(includes CCHO)FF 




chalcomycin- 
myconomycin- 
aldgamycinDmiko 
nomycm 


BNNGKDN 


B(2-0-glycosyl)N(l,2-epoxy)NG(2- 
hydroxyl)KD( 1 -glycosyl)N 


S. bikiniensis, S. 
albogriseolus 


chimeramycinB- 
FTL448 


BMNCODF 


B(0-ethyl)MNC(l -glycosyl)OD(l - 
glycosyl)F 


S. ambofaciens ka- 
448 


chivosazolA 


SRRnNNSRQNn 
RMnNn 


SRR(l-0-macrocyc)nNNSR(l- 

methoxy)QNnR(l- 

glycosyl)MnNnQ(keto-macrocyc) 


S. cellulosum 


cineromycinB 


BLME=JN 


BLME=J(2-hydroxyl)N 


S. 

cinereochrornogene 
s, S. sp. 


cineromycinBdehy 
ro 


BLMI=JN 


BLMI=J(2-hydroxyl)N 


S. grieoviridis 


cineromycinB2,3di 
hyro 


BLME=JL 


BLME=J(2-hydroxyl)L 


S. grieoviridis 


cirramycinB 1 dihy 
droxy-A6888X 


BMNGODF 


BM(l > 2-epoxy)NGOD(l-glycosyl)F 


S. flocculus 


cirramycinB - 
cirramycinB 1- 
Acumycin- 
A688A-B58941- 
A6888A 


AMNGODF 


AM(l,2-epoxy)NGOD(l -glycosyl)F 


S. cirratius, S. 
griseoflavus, S. 
fradiae, S. 
flocculus 


cladospolideA 


ELLFN 


ELLF(2-hydroxyl)N 


Cladosporium 
fulvum, C. 
cladosporiodes 


cladospolideB 


ELLFn 


ELLF(2-hydroxyl)n 


Cladosporium 
fulvurn, C. 
cladosporiodes 


cladospolideC 


ELLEN 


ELLE(2-hydroxyl)N 


fungus 

Cladosporium 
tenuissimum 
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POLYKETIDE 


CHUCKLES 


annotated-CHUCKLES 


SOURCE 


concanamycrnA- 
folimycin-A661 -1 - 
S45A-TAN1323B- 
X4357B 


CENMKCCMN 


CE(2-methoxy)NMKC(2- 
ethyl)CMN(2-methoxy) 


S 

diastatochromogen 
ese, o. sp, o. 
neyagawaensis 


concanamycinB- 
S45B 


CENMKCCMN 


CE(2-methoxy)]^MKCCMN(2~ 
methoxy) 


s 

diastatochromogen 
ese 


concanamycinG- 

anhydroconcanam 

ycinB 


NBNHCENMKC 
CMN 


NBNHCE(2- 

rnethoxy)NMKCCMN(2-methoxy) 




congloblatin 


AJM 


A(l-oxazoyl)JM.AJM 




#copiamycin 


RNRSSSSSSRLR 
SRN 


RNRSSSS(l-0-cyc)S(2- 
hydroxyl)S(l-cyc)RLRSRN 


S. hygroscopicus 


cytovaricin-H230 


QQQKQNLJFCD 
N 


QQQKQNLJ(2-hydroxyl)F(2-0- 
glycosyl)CD(2-hydroxyl)N 


S. sp., S. collmus 


cytovaricinB 


QQKQNLJFDDN 


QQQKQNLJ(2-hydroxyl)F(2-0- 
glycosyl)CD(2-hydroxyl)N 


S. torulosus 


CP64537 


ARGKDD 


A(2-glycosyl)R(l- 
C(=0)C(0)C(C)QG(2- 
hydroxyI)KDD( 1 -glycosyl) 


Streptomyces 
toyocaensis 
humicola A ICC 
39491 


damavaricinC 


QDQCCNM 


QDQCCNM 


S. spectabilis 


deltamycinAl 


ENNGOFF 


EN( 1 ,2-epoxy)NGOF( 1 -giycosyl,2- 
methoxy)F(l-C(=0)C) 


a. deltae, a. 
halstedii-deltae 


deltamycinX- 

desisovalerylcarbo 

mycinA 


ENNGOFF 


EN(l,2-epoxy)NGOF(I-glycosyl,2- 
methoxy)F(l-C(=0)C) 


S. deltae, S. 
halstedii-deltae 


engleromycin 


QNJHN 


QNJH(2-hydroxyl)N(U-epoxy) 


Engleromyces 
goetzei 


#epothilone 


MEMLJDgE 


MEM(1 ,2-epoxy)LJDG(2~rnethyl)E 




erythromycin 


ADGJDD 


A(2-hydroxyl)DGJ(2-hydroxyl)D(l-glycosyl)D(l- 
glycosyl) 


espinomycinA2 


ENNCOFF 


ENNCOF(l -glycosyl,2-methoxy)F(l - 
C(=0)C) 


S. fungicidicus 


filipinlH- 
lagosinl4deoxy 


ENNNNMEFFFF 
FFFL 


E(2%droxyl)NNNlSnVfEFFFFFFF 


S. fllipinensis, S. 
durhamensis 


filipin-lagosin 


FFF 




S. fllipinensis, S. 
durhamensis 


formamicin 


CBNMOCMN 


CBNMO(includes a long, branched 
alkyl chain)CMN 




foromacidinB- 

spiramycinE- 

spiramycinB 


ENNAOFF 




S. ambofaciens 


FD891 


RRUSNNLNRM 
M 


RRUS(1 -0-macrocyc)NNL(2- 
hydroxyl)N(l ,2-epoxy)RMMQ(keto- 
macrocyc) 


S. grarninofaciens 


FK895 


RNRNMRNRLS 


R(l -methoxy)N( 1 ,2-epoxy)RNMR(l - 
0-rnacrocyc)NR(l -C(=0)C,2- 
hydroxyl)LSQfketo-macrocyc) 


S. hygroscopicus 


FK-506 


MAEPMJJBKQQ 






gedamycin 


IEJBNNNnNNNF 
AEFEEEEIE 


IEJBNNNriNNNF(l»glycosyl)A(l-0- 

cyc,2-carboxylicacid)EF(2- 

cyc)EEEEIE 


S. aureofaciens 
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PfYI WTTTTTIT? 
i UL x iVJL 1 11/ JL 






ORGANISM 


geldanamycin 


QULRMRnM 


QUL(2-metnoxy)RMS(l -CONH2,2- 
methoxy)nM 


S. hygroscopicus 
var. gclanus 


gephyronicacid 


RRTSRRM 






gloeosporone 


ELLEEL 


ELLE(l-0-cyc)I(2^cyc^-hydroxyl)L 


Colletotrichum 
gloeosporioides f. 
sp. jussiaea 


GERI-155 


BNLGJDN 


B(2-0-glycosyl)N(l,2^poxy)LG(2- 
hydroxyl)JD(l -glycosyI)N 


S. GERI-155 


halomicin 


QNCRCBNM 


QNC(1 -methoxy)R(l - 
C(=0)C)CBNM 


Mic. halophytica 


herbimycin 


ALDMFnM 


A(l-methoxy)L(2-methoxy)D(l- 
methoxy)MF(l -CONH2,2- 
methoxy)nM 




hygrolidin 


CENMJCMM 


CE(l-0-macrocyc,2- 

methoxy)NMJCMMQ(keto- 

macrocyc) 


S. hygroscopicus 


hygrolidin-oxo 


BNHCENMJCM 
N 


BNHCE( 1 -0-macrocyc,2- 

methoxy)NMJCMN(2- 

methoxy)Q(keto-macrocyc) 


S. griseus, S. 
hygroscopicus 


irumamycin 


URRNSLMQQS 


UR( 1 -0-macrocyc)RNS( 1 - 

glycosyl)LMR(l-0-cyc)=QS(2- 

cyc)Q(keto-macrocyc) 




juvenimicinAl- 

T1124A1- 

M4365A1 


BMNGFDF 


BMNGFDF(includes an ethyl in pos 
2) 


Mic. chalcea 


juvenimicinA2- 

T1124A2- 

M4365A2 


BMNGJDF 




Mic. chalcea 


juveniniicinA4- 
Tl 124A4- 
M4365A4 


BMNGODF 




Mic. chalcea-Mic. 
capillata 


juvenmucin-T1124 


BNNGJDF 




Mic. chalcea 


kanchanamycin 


NMSSSSSSSRL 
RSRNN 






lankamycin- 

kujimycin- 

landavamycm- 

a ono O OX TO 

A20338N2 


ADGJDD 




S. violaceoniger, S. 
spinichromogenes 


leinamycin 


QNN1MLN 






leucanicidin 


CENMKCMN 


CE(2-methoxy)NMKCMN(2- 
methoxy) 


S. halstedii . 


leucomycinA12- 
kitasomycinA12 


FNNCODF 




S. kitasatoensis 


leucomycinA14- 
kitasomycinA14 


FnNCOFF 




S. kitasatoensis 


leucomycinA3- 
josamycin- 
platenomycinA3 - 
turimycinA5 


!■' V TXT A /"M^vT** 

ENNAODF 




S. kitasatoensis, S. 
hydroscopicius, S. 
narbonensis, S. 
platensis 


leucomycinAS- 
taimcinH4 


ENNCOFF 




S. kitasatoensis 


lienomycin 


SSSNSSSTSML 
NRRNNNNNL 
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POLYKETIDE 


CHUCKLES 


annotated-CHUCKLES 


SOURCE 
ORGx\NISM 


lucensomycin- 
etruscomy cin - 
lucirnycin-FJ 1 1 63- 
butylpimaricin 


SNNNNSREEFN 
N 




Act. sp, S. lucensis, 
S. glaucus 


L155175 


RSNMURmM 


RS(2-methoxy)NMURmM 




L681110 


NMURMN 


NMURMN(2-methoxy) 




macrocin- 
lactenocin 


BMNHODF 




S. aureus, S. lutea, 
K. pneum, B. 
subtl., Shva.. 


macrocin-Y07625 


QMNGODF 




S. fradiae gs 16 


maridomycinl 


ENNCOFF 




S. paltensis, S. 
ritnosius, S. 

<"*afn)PTi<!iQ R 

racemochromogene 

s 


maridomycm- 
platenomycinC3- 
ti irinrvcinFP S - 

LUl 1111 jvIIJIjI J 

B5050A-YL704- 
C-3 


ENNCJDF 




S. hygroscopicus- 
S. platensis- 
malvinus 

XXXU.J » XXX uo 


TnatfiPTnvrin A 


RSSRSSRSSRR 
MRMRNLRL 






midecamycinAl- 
platenomycinBl- 
SF837 


ENNCOFF 


ENNCOF(2-metboxy)F 


S. mycarofaciens 


mideeamveinA2- 
mydecamycinA2- 
SF837A2 


FNNDJEE 


FNND JEf2-hvdroxv^Ef 1 -C(=U)CO 


S mvcarofaciens 


milbemycin 


QQQMKNQQ 






#monazomycm 


SSRRSSURNSR 
SMRMRNLRSL 
L 






m vein ami cin 

x 1 xjr wn xcxxx xx w ju i 


BNNGJDN 


B ( 1 -c vcYNNG JDN(2-cvcl 




mvcin ami rin VT 


BNNGJDN 




Micromonosnota . 
griseorubida sp. 


mycinamicinXl 1 


BNNGLDN 




S. aureus, S. 

ovoeenes, 

Corvnebacterium 


#my c olactone 


SRMUSMSSL.sS 
MNMMN 


SRMUSMSSL.sSMNMMN 




mycolactoneA 


SRMUSMSSL 


SRMUSMSSL 




mycolactoneB 


sSMNMMN 


sSMNMMN 




mvxovire^cinA 1 

XXX Y **.\J T J A WOWllU V X 


OOEFLNNLL IT , 
KJ 






mvxovirpscm A 9 


OOEFOJNLLILJ 

J 






myxovirescinJB- 
megovalicinB 


QQEFLNNLLIL 
KM 






myxovirescinC-C 1 


QQEFLNNLLLL 
KJ 






myxovirescinD 


QQEFLNNLLLL 
KM 






myxovirescinE 


QQEFLNNLLDLJ 
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POLYKETIDE 


CHUCKLES 


annotated-CHUCKLES 


ISUUKCll/ 
ORGANISM 


myxovirescinFl 


QQEFLNNLLLL 
KJ 






myxovirescinF2 


QQEFLNNLLLL 
JJ 






myxovirescinGl 


QQEFLNNLLLK 
J 






myxovirescinG2 


QQEFLNNLLLJJ 






myxovirescinHl 


QQEFLNNLLIL 
KJ 






myxo vire scinH2 


QQEFLNNLLILJ 
J 






myxovirescinL 


QQEFLNNLLIL 
HJ 






myxovirescinP 1 


QQEFLNNLLIL 
LKJ 






myxovirescinP2 


QQEFLNNLLIL 
LJJ 






myxoviiescinQ 


QQEFLNNLLIL 
KM 






myxovirescinS 


QQEFLNNLLIL 
HJ 






M4365G2 


BMNGODF 




Streptoverticillum 
kitasatoensis, S. 
thenhaotolerans 


nancimycin 


QNDACBNM 




S. albovinaceus 


neocopiainycin 


NRSSSSSRLRSR 
N 






niddamycin- 
F3463- 

3 desacetyjcarbomy 
cinB 


FNNGOFF 


FNNGOF(l -glycosyl,2-methoxy)F 


S. aureus, S. lutea, 
B. subt 


oligomycinA 


QJNNJRGAGAN 




diastatochromogen 
es, S. chibaensis 


oligomycinB 


QJNNKCHDHD 
N 




S. 

diastatochromogen 
es 


oHgomycinB- 
44homo 


QJNNJRGRGAN 




S. bottropensis 


oligomycinD 


QJNNJBGAGAN 




S. arabicus, S. 
parvulus, S. 
rutgersensis, S. 
griseus, S. 
aureofaciens 


ossamycin 


QQQNLLKFDD 
N 






perimycin 


JBNNNnNNNFC 
EFEEEEIE 


JBNNNnNNNFCEFEEEEEE 




phenalamid 


nCNNNNM 






phenalamideAl- 
fenalamid-I-102-C 


JMCMNnNM 






phenal amide A2- 
102-T 


JMCMNNNM 






phenalamideA3 


JMCMNNnM 







WO 01/92991 



28 



PCT/US01/17352 



POLYKETIDE 


CHUCKLES 


annotated-CHUCKLES 


SOURCE 

UJtVvjr/\i> lolVI 


. — 

phenalamideB 


JMCJMJNnJNM 






pnenalamiaeC 


TN A ft ifXTXTKTN >f 

JMCMJNINjnJVI 






phthoramycin 


r\fwn T TOD "D XT 






pikromycin 


ANGJDH 






prasinon-L155175 


DNMJCMM 




S. hydroscopicus 
ma jZod, 0. 
prasinus 


protostreptovaricin 


QDACDNM 






PD118576A2 


T7VT* JTOI XTlT 

ENMJCMN 




0. sp. wp 3913 


PF1163 


EJLLQ 






tfquinolidomicin 


SSLULSSNNNSI 
SNNSSSUSRRQ 
QRRSS 






rapamycin 


FGMEGJNNME 
EKQQ 






rhizopodin 


RLQSNNSS 






rifamycin 


QQQNDACBNM 


— , 1.,,,. v 




rimocidin 


BNNNNFCEFEL 
IA 






rosamicin- 
repiomicin 


BMNGODF 






rosaramycin 


BMNGODF 




Micromonospora 
rosaria 


rustmicin 


QQQJMOG 


QQQJMCKincIudes COH)G 




rutamycin 


QQJNNJBGAGA 
N 






scytophycin 


BFCEENQEMN 






scytophycinB-E 


BFCEQNQEMN 






shurimycin 


NNRSSSSSSRLR 
SRNN 






sorangicin 


LKMFnENLNCF 
DFNQNNnn 


LKMF(2- 

hydroxyl)nENLNCFDFNQNNiin 




sorangicinA 


QNLNF 






sorangicinB 


NLNF 






sorangolideA 


LLLLKBUMEA 
FM 




myxobacterium 

sorangium 

cellulosum 


soraphen 


TIT TpfTN'p A 

ELJFJDFA 






spiramycin 


nNCJDF 






staphcoccomycin- 

angolamycin- 

shincomycin 


CMNGODF 






stipiamide 


JMCNNnNM 






tartrolonBz 


XTT tTTT? 

nNLFHE 


nNLFHE 


— — ; 

fragment 


tedanolide 


JGEHMDHF 


JGEHMDHF(2-hydroxyl) 




thiazinotrienomyci 
n 


r"\T _~T) CTsTNTXTC 






tiacumicin 


SMMRMSNM 






tylosinC-macrocin 


EMNHODF 






tylosin-A 


EMNGODF 






TAN-1323 


NMUREMM 






TMC-34 


NRSSSSSSRLRS 
RN 
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POLYKETIDE 


CHUCKLES 


annotated-CHUCKLES 


SOURCE 


venturicidinA 


CKNFLMQQQ 






vicenistatin 


LNMULRNN 






viranamycinB- 

virustomycin- 

TAN1323C 


QNMKCCMN 




S. sp. ch41 


zmcophorin- 
griseochelin 


KMCNLACC 




S. griseus 



In another embodiment, the polyketide library includes the name of the polyketide, the 
CHUCKLES string and a linearized representation of the structure. The linearized 
representations of the CHUCKLES structures for erythromycin and epothilone are as follows: 



epothilone D 
MEMUDqE 



erythromycin 
ADGJDD 




An illustrative example of a polyketide library containing linearized representations of their 
structures is found in Appendix C (deposited in the microfiche appendix). 

EXAMPLE 3 

Alternative PKS Genes for Epothilone 

This example illustrates the alignment and design of novel PKS genes for the target 
epothilone. Epothilone is first converted into CHUCKLES string format and then read into the 
MORPH program as a TARGET. The program then generates all possible alignments of library 
modules and sorts the alignments to determine preferred combinations of modules for gene 
construction and production of epothilone via a novel polyketide synthase gene. 
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O OH 0 



The epothilonc D structure above was first opened at the macrolactone ring closure 
between the C-l-ketone and the C-15-oxygen. The monomer set shown in Figure 2 was then 
matched against each of the successive pairs of macrocyclic backbone carbon atoms, starting 
with C-2 and C-3, which match monomer E. The next two carbon atoms C-4 and C-5 match 
monomer G with an additional post-synthetic methylation on C-4. C-6 and C-7 match monomer 
D. C-8 and C-9 match monomer J. C-10 and C-l 1 match monomer L. C-12 and C-13 match 
monomer M. C-14 and C-l 5, where C-15 has a hydroxyl substitution (modified by thioesterase 
to close the macrocycle), match monomer E. C-l 6 and C-l 7 match monomer M. 

The rest of the molecule, a methyl-substituted thiazole moiety, does not match any of the 
monomers in the monomer set. This moiety corresponds to a malonyl CoA loading module and 
an NRPS module that together generate the methyl-substituted thiazole moiety. This moiety is 
thus omitted from the CHUCKLES string generated from this illustrative monomer set but can 
be added simply by adding a monomer to the set. The CHUCKLES string generated is 
EGDJLMEM, which is in the reverse order of biosynthesis. This sequence is then reversed to 
MEMLJDGE to yield a monomer sequence that matches the order of biosynthesis. The 
sequence is then annotated to account for the post-synthetic modifications as follows 
MEMLJDG (2-methyl)E. 

This target sequence is provided to the MORPH program to generate all possible 
combinations of modules in the CHUCKLES-encoded library that will yield the target 
CHUCKLES. The valid combinations are then sorted in increasing order of non-native inter- 
module interfaces. In one implementation, a MORPH run generated 3,452 valid sequences of 
five inter-module interfaces. Of these, none contain fewer than five inter-module interfaces. 
Some illustrative sample module combinations appear below. The combinations are shown 
listing each monomer followed by a colon and the name of the polyketide(s) from which it is 



WO 01/92991 



31 



PCT/US01/17352 



derived, followed by a parenthetical showing the associated monomers, in that polyketide. 
Vertical lines represent modular junctions between two different polyketides. 

Illustrative PKS Gene 1: 

MrSacetyW'TjutylryltylosinfFMN) | Ertedanolide(GEH) | M:aldgamycin(BML) 
L:aldgamycin(MLG) | J:aldgamycin(GJD)D:aldgamycin(JDL) | G:tedanolide(JGE) 
E:tedanolide(GEH) 

Illustrative PKS Gene 1 thus comprises one or more open reading frames that encode, in 
the order listed, the module from the acetyI~4"-butyryltylosin PKS that corresponds to monomer 
M, the module from the tedeanolide PKS corresponding to monomer E, the modules from the 
aldgamycin PKS corresponding to monomers M, L, J, and D, and the modules from the 
tedanolide PKS corresponding to monomers G and E. 

Illustrative PKS Gene 2: 

M:albocycline-Ml-ingramycin-TA2407-cineromycinB-U28010-SR20 
Ml- ingraniycin-TA2407-cineromycinB-U28010-SR2077(MEJ) | M:albocycline-Ml- . 
ingramycin-TA2407-cineromycinB-U28010-SR2077(LME) | L:albocycline-Ml- ingramycin- 
TA2407-cineromycinB-U28010-SR2077 (BLM) | J:erythromycin(GJt)) D:erythromycin(JDD) | 
G:tedanolide(JGE) Ertedanolide(GEH) 

EXAMPLE 4 

Alternative PKS Genes for 6-Deoxvervthronolide B 

This example illustrates the alignment and design of novel PKS genes for the 
erythromycin basic polyketide structure (6-dEB) using the MORPH program. 
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For the 6-dEB structure above, the CHUCKLES string is generated by first opening the 
macrolactone ring closure between the C-l-ketone and the C-13- oxygen. Using the monomer 
set and matching protocol described in Example 3, one generates the CHUCKLES string 
DDJGDA, in the reverse order of biosynthesis. This sequence is then reversed to ADGJDD to 
yield the monomer sequence that matches the order of biosynthesis. The sequence is then 
annotated to account for the post-synthetic modifications (erythromycin A) as follows A(Z- 
hydroxyl)DGJ(2-hydroxyl)D(l-glycosyl)D(l-glycosyl). 

This target sequence is supplied to the MORPH program to generate all possible 
combinations of modules in the CHUCKLES-encoded library. The valid combinations are then 
sorted in increasing order of non-native inter-module interfaces. In one implementation, a 
MORPH run generated 19,631 valid sequences of less than or equal to five inter-module 
interfaces. Of these, 13,306 contain 4 inter-module interfaces, and 256 contain only 3 inter- 
module interfaces. Some of these contain only two inter-module faces, and one only contains 
one. Some illustrative sample module combinations follow. 

Illustrative PKS Gene 1 : 

A:amphotericinA(EAL) j D:aldgamycin(JDL) | G:mycinamicin(NGJ) J:mycinamicin(GJD) 
D:mycinamicin(JDN) | D:amphotericinA(CDN) 

Illustrative PKS gene 1 thus comprises one or more open reading frames that encode, in 
the order listed, the amphotericin PKS module corresponding to monomer A, the aldgamycin 
PKS monomer corresponding to monomer D, die mycinamicin PKS modules corresponding to 
monomers G, J, and D, and the amphotericin PKS module corresponding to monomer D. 
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Illustrative PKS Gene 2: 

A:amphotericinA(EAL) | D:aldgamycin(JDL) | G:pikromycin(NGJ) J:pikromycin(GJD) 
D:pikromycin(JDH) | D:aldgamycin(JDL) 

Illustrative PKS Gene3: 

A:larikamycin-kujiinycin-landavamycin-A20338N2(-AD)Dd 
landavamycin-A20338N2(A^ 

J:lankamycin~kujimycin-landavamycin-A20338N2(GJD) | D:ossamycin(FDD) 
D:ossamycin(DDN) 

Illustrative PKS Gene 4: 

A:amphotericinA(EAL) | D:lankamycin-laijimycin-landavamycin-A20338N2(ADG) 
G:lankamyck-kujimycin-landavamycin-A20338N2(DGJ) 
landavamycin-A20338N2(GJD) | D:A82548A-cytovaririn(FDD) D:A82548A- 
eytovaricin(DDN) 

Illustrative PKS Gene 5: 

A:lankainycin-kujirnycm-landavamycin~A20338N2(-AD) D:lankamycin-kujimycin- 
landavamycin-A20338N2(ADG) G:larjkamycin-kujimycin4andavamycin-A20338N2(DGJ) 
J:lankarnycin-kujimycm-landavamycm-A20338N2(GJD) Dilankamycin-kujimycin- 
landavamycin-A20338N2(JDD) D:lankam^ 

EXAMPLES 

Source Code: 
#include <stdio.h> 

/* -siani/programs/morph/morph3.c 

PURPOSE: To traverse recursively all the entries in PKS.lib, 
generating all feasible combinations of PKS modules to make the TARGET (e.g., epothilone). 
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INPUT: libraryfile: tab-delimited CHUCKLES-coded polyketides file with the 
following columns: 

1 . polyketide name 

2. plain CHUCKLES 

3. annotated CHUCKLES (contains information about post- 
synthetic modifications) 

4. source organism; 
targetname: user-defined name (e.g., epoD); 

targetsequence: CHUCKLES-coded polyketide of desired TARGET (e.g., 

MEMLJDGE); 

X, Y, Z sets of wildcards: sets of monomers for particular positions 
appearing in target sequence (the wildcards can be used for analoging the TARGET polyketide); 
hard-coded parameters which may be reset (requires recompiling): 

NBOUNDARY_CUTOFF determines the maximum number of 
non-native inter-modular interfaces which are contained in the output (set to 5, but may be 
increased or decreased); and 

RECURSION J30UNTER_CUTOFF specifies the number of 
levels of recursion (defaults to 0, 1, 2) acceptable for the run — a large PKS library can cause 
recursion that will greatly increase run time; because of the multi-directionality of the 
alignments (using every library entry as a STARTER), there is typically no need to go beyond 2 
levels of recursion. 
OUTPUT: 

All combinations of modules that meet parameters set by user. Example 
output from MEMLJDGE (epothilone D) using a subset of a PKS library is provided below. 
Vertical bars indicate non-native inter-modular interfaces. Last column contains the number of 
"pieces" that are needed to put together the PKS. 

Names of PKSs have been abbreviated to fit them in these 

comments. 

HIT M:3atyl(FMN)| E:tedan(GEH)| M:aldga(BML) L:aIdga(MLG)| 

J:aldga(GJD) D:aldga(JDL)| G:tedan(JGE) E:tedan(GEH)| 5 

HIT M:albMl(LME) E:albMl(MEJ)| M:albMl(LME)| L:aldga(MLG)| 

J:aldga(GJD) D:aldga(JDL)| G:tedan(JGE) E:tedan(GEH)| 5 

HIT M:albMl(LME) E:albMl(MEJ)| M:aldga(BML) L:aldga(MLG)| 

J:aldga(GJD) D:aldga(JDL)| G:3aty](NGO)| E:albMl(MEJ)| 5 

HIT M:albMl(LME) E:albMl(MEJ)| M:aldga(BML) L:aldga(MLG)| 
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J:aldga(GJD) D:aldga(JDL)| G:aldga(LGJ)| E:albMl(MEJ)| 5 
HIT M:albMl(LME) E:albMl(MEJ)| M:aldga(BML) L:aldga(MLG)| 
J:aldga(GJD) D:aldga(JDL)| G:aldga(LGJ)| E:albMl(MEJ)| 5 
USAGE: 

morph3 -1 libraryfile -n targetname -t targetsequence [-x X-wildcards] [-y Y- 
wildcards] [-z Z-wildcards] 
examples: 

# generate combinations that yield epothilone D 

%morph3 -1 PKS.lib -n epoD -t MEMLFDGE > omorph3_epoD 

%egrep HIT omorph3_epoD | sort | uniq | sort +10 -1 1 > 
omorph3epoD.uniq.sort 

%egrep ALIGN JTARGET omorph3_epoD > 
omorph3_epoD_STARTER_ALIGN 

. # generate combinations that yield epothilone D with a 
C13-hydroxyl 

%morph3 -1 PKS.lib -n epoD-130H -t MEXLJDGE -x ABCD > 

oepoD-130H 

%egrep HIT oepoD-130H | sort | uniq | sort +10 -1 1 > 
oepoD- 1 3 OH . uniq .sort 

%egrep ALIGN_TARGET oepoD-130H > oepoD-130HJ5TARTER_ALIGN 

# generate combination that yield epothilone with the 
following wildcards (set 1) 

%morph3 -1 PKS.lib -n epoD-setl -t MEXYZDgE -x ABCD -y LEFIN -z 
JACGM > oepoD-setl 

%grep HIT oepoD-setl | sort | uniq | sort +10 -1 1 > oepoD-setl. uniq. sort 

# generate combination that yield epothilone with the following wildcards (set 2) 
%morph3 -1 PKS.lib -n epoD-set2 -t MEXYZDgE -x JK -y EF -z JACGM > 

oepoD-set2 

%grep HIT oepoD-set2 ) sort | uniq | sort +10 -1 1 > oepoD-set2.uniq.sort 
LIMITATIONS: 

This version does not handle intra-modular modifications/splitting because 
morph is operating at the monomer level. Modifications could convert the CHUCKLES-encoded 
strings into the corresponding and equivalent SMILES and then perform more complex 
chemical analysis of the PKS molecular graphs. 
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Currently, inter-modular double bonds are present in the library, but are ingnored by the 
morph program. 
*/ 

#include <stdio.h> 

/* ~siani/programs/morph/morph3.c 
*/ 

#defineTRUE 1 
#defineFALSE 0 
#defineDEBUG_MATCH FALSE 
#defmeDEBUG_STARTER FALSE 
#defineDEBUG_ALIGN FALSE 
#defineDEBUG_RECURSE FALSE 
#defineDEBUG_WILDCARD FALSE 
#defineMAXLEN 80 
#defineMAX_TYL_LEN 6 
#defineMAX_EPO_LEN 6 
#defineMAXNAMELEN 160 
#defineMAX_LIB_ENTRIES 500 
#defineMAXWTLD 3 
#defineMAXBUF 200 
#defineNB OUND AR YCUTOFF 5 
#defineRECURSION_COUNTER_CUTOFF 2 
#defineSTARTER_MINIMUM_ADJACENT_ALIGN 2 
MefinelVnNIMUM_ADJACENT_ALIGN 2 
typedef struct _lib { 

char name[MAXNAMELEN]; 

char monomersequence[MAXNAMELEN] ; 

char annotatedsequencefMAXNAMELEN]; 

char aUgnedsequencefMAXNAMELEN]; 

char alignedPKSname[MAXLEN] [MAXNAMELEN]; 

int boundarytoright[MAXNAMELEN]; 

int marked[MAXLEN]; 

char context[MAXLEN][4] ; 

int recursiontagged; 

int nboundary; 
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}LIB; 

main(int argc, char **argv) 
{ 

int ii=0 Jj=0, kk=0, 11=0; 

int nlib=0; 

int ecounM); 

int nfilled=0,nfilledmax=0; 

int epothilonelen=0; 

int nlargestpiece==0,tnlargestpiece=0; 

int mrnpass=0; 

int lcount=0; 

int new_unmarked_entries_filled=0; 
int recursion_counter = 0; 
int nwildcard=0; 

int best_new_unmarked entries filled = 0; 

int smallest_acceptable_j>iece = 0; 

int current_nmarked=0,previous_nmarked^O; 

char *sptr, *eptr, *Iptr,*bufptr; 

char *clibptr; 

char *libraryfile; 

char *targetsequence,*targetname; 

char bui[MAXBUF]; 

char wildcards[MAXWIIJ)]fMAXLEN]; 

FILE *libp; 

LIB epotemp; 

LIB library[MAX__LIB_ENTRIES]; 

LIB epothilone; 

char *progname; 

char **filelist, **fileptr; 

libraryfile = f,H ; 

targetsequence = n "; 

targetname= 

for(ii=0; ii<MAXWILD; ii++) { 

for(jj=0;jj<MAXLEN; jj++) { 
wildcards[ii][jj] = '\0'; 
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} 

} 

/* process arguments */ 

filelist = fileptr = (char **)(malloc(argc * sizeof(*argv))); 
progname = *argv++; 
if(argc<2) { 

§)rintf(stderr, ,f usage:%s -1 libraryfile -n targetname -t targetsequence [~x X- 
wildcards] [-y Y-wildcards] [-z Z-wildcards]\n M ,progname); 
exitO; 

} 

while(argc~ > 1) { 

if(argv[0][0] == && argv[0][l] != \W) { 
/* handle option */ 

*-H-(*argv); /* advance past the minus */ 
switch(**argv) { 

case T: /* get library input filename (PKS.lib) */ 

argv-H-; argc-s 

libraryfile = argv[0]; 

fprintf(stderr, rt -l: libraryfile=%s\n^libraryflle); 
break; 

case V: /* get target name string */ 
argv-H-; argc--; 
targetname = argv[0]; 

fprintf(stderr, w -t: targetname=%s\n",targetname); 
break; 

case V: /* get target sequence string */ 
argv-H-; argc-; 
targetsequence = argv[0]; 

fprintf(stderr,"-t: targetsequence=%s\n , \targetsequence); 
break; 

case V: /* get a wildcard string */ 
argv-H-; argc-; 
strcpy(wildcards[0],argv[O]); 
fprintf(stderr,"-x: 
wildcards[%d]=%s\n n ,0,wildcards[0]); 
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nwildcard++; 
break; 

case y : /* get a wildcard string */ 
argv-H-; argc--; 
strcpy(wildcards[l],argv[0]); 
fprintJfitstderr/'-x: 
wildcards[%d]=%s\n M > l,wildcards[l]); 

nwildcard++; 
break; 

case % z': /* get a wildcard string */ 
argv-H-; argc-; 
strcpy(wildcards[2],argv[0]); 
fprintf(stderr, M -x: 
wildcards[%d]=%s\n n ,2,wildcards[2]); 

nwildcard-H~; 
break; 

case 'X 1 : /* get a wildcard string */ 
argv-H-; argc-; 
strcpy(wildcards[0],argv[0]); 
fprintf(stderr,"-x: 
wildcards[%d]=%s\n ,, ,0 5 wildcards[0]); 

nwildcard-H-; 
break; 

case 'Y': /* get a wildcard string */ 
argv-H-; argc-; 
strcpy(wildcards[ 1 ],argv[0]); 
fpriQtf(stderr,"-x: 
wildcards[%d]=%s\n ,, ,l ,wildcards[l]); 

nwildcard++; 
break; 

case r Z': /* get a wildcard string */ 
argv-H-; argc-; 
strcpy(wildcards[2],argv[0]); 
fprmtf(stderr, n -x: 
wildcards[%d]=%s\n ,, ,2,wildcards[2]); 
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nwildcard-H-; 
break; 
default: 

fprintf(stderr,"%s unknown option; 

ignored\n'\*argv); 

}/*switch*/ 
} else { /* a regular filename */ 
*fileptr++ = *argv; 
*fileptr = NULL; 

} 

argv-H-; 

}/*while*/ 
if(nwildcard > 0) { 

for(ii=0; ii< nwildcard; ii-H-) { 
^rintf(stderr/wUdcards[%d]=%s\n ,, J ii J wildcards[ii]); 
fyrintf(stdout;Vildcards[%d]=%s\n ,, ,ii > wildcards[ii]); 

} 

} 

epothilone Jiboundary = 0; 
for(ii=0; ii<MAXNAMELEN; ii++) { 
epothilone.name[ii] = 'VO'; 
epothilone.monomersequence[ii] = '\0'; 
epothilone.alignedsequence[ii] = W; 
epothilone.boundarytoright[ii] = TRUE; 
forGj=0; jj<MAXLEN; jj++) { 

epothilone.alignedPKSname[jj][ii] = \0 U , 

epothilone.marked[jj] = FALSE; 

epothilone.context[j[j][0] = '\0 ! ; 

epothilone.context(jj][l] = f \0'; 

epothilone.context[jj][2] = 

} 

} 

strcpy(epothilonejiame,targetname); 
strcpy(epothiloncjnonomersequence,targctsequence); 
fprintf(stdout, "TARGET: %s\n", epothilone.monomersequence); 
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ecount = 0; 

eptr = epothilone.monomersequence; 
while(*eptr!-W){ 

if(ecount = 0) { 

epothi]one.context[ecount][0] - ; 

epothilonexontext[ecount][l] = *eptr; 

epothilone.context[ecount][2] = *(eptr + 1); 

} else { 

if(ecount = 
(strlen(epothilone.monomersequence) - 1)){ 

epothilone.context[ecount][0] = 



*(eptr-l); 



*(eptr-l); 



*(eptr+l); 



epothilone.context[ecount][l] = *eptr; 
epothilone.(X)ntext[ecount][2] = ; 

} else { 

epothUone.context[ecount][0] = 

epothilone.context[ecount][l] = *eptr; 
epothilonc.contexttecoxint][2] = 



} 

} 

epothilone.context[ecount][3] = W; 

eptrf+; 

ecount++; 

} 

for(ii=0; ii<ecount; ii++) { 

^>rintf(stdout,"(%s)\n",epothilone.context[ii]); 

} 

/* library ■ */ 

nlib = get_library(libraryfile,library); 

fyrintf(stdout,"nlib=%d\n ,, ,nJib); 

kk=0;while(kk<nlib){ 

/* zero out the epothilone entry with respect to a new alignment • 

for(ii=0; ii<MAXNAMELEN; { 
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epothilone.aiignedsequencefii] = \Q'\ 
epothilone.boundarytoright[ii] = TRUE; 
for(jj=0; jj<MAXLEN; jj-H-) { 

epothilone.alignedPKSname[jj][ii] = '\0 f ; 

epothilone.marked[jj] = FALSE; 

} 

} 

/* - reset the context back to that in epothilone- */ 
ecount = 0; 

eptr = epothilone.monomersequence; 
while^eptrN'VO'K 

if(ecount = 0) { 

epothilone.context[ecount][0] = ; 

epothilone.context[ecount][l] = *eptr; 

epothilone.context[ecount][2] = *(eptr + 1); 

} else { 

if(ecount = 
(strlen(epothilone.monomersequence) - 1)){ 

epothilone.context[ecount][0] = 

*(eptr-l); 

epothilone.context[ecount][l] = *eptr; 
epothilone.context[ecount][2] = 

} else { 

epothilone.context[ecount][0] = 

*(eptr-l); 

epothilonexontext[ecount][l] = *eptr; 
epothilone.context[ecount][2] = 

*(eptr+l); 

} 

} 

epothilone.context[ecount][3] = \0 tm 9 

eptrH-; 

ecount*H-; 

} 

/* align STARTER (current library entry) and 
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epothilone */ 

sptr = libraryfkk] .monomersequence; 
lcount=0; 

while(*sptr != W) { 

^rintf(stdout/1ibrary[%d].monomersequence[%d]==%c\n", 
kkJcount,Ubraiy[kk].monomersequence[Icoiint]); 

sptr++; 
lcount-H-; 

/* Call maximal_adjacent_alignment until it no longer 
returns more than two adjacent modules. There is no reason to 

try to extract individual modules, because this is done as part of the recursive filling of spaces 
from the library. 

*/ 

smallest_acceptable_piece = 2; 
eptr = epothilone.monomersequence; 
flprintftstdout, " ALIGN JTARGET: "); 
while(*eptr != 'VO') { 

fprintf(stdout,"%c" J *eptr); 

eptrH-; 

} 

$rintf(stdout/ f \n"); 

fprintf{stden-, M aligning %d %s\n",kk, library[kk].name); 
best_new_unmarked_entries_filled = 0; 
while((new_unmarked_entries_filled = 
maximaladjacentalignmentanddmn^ 

mallest_acceptable_piece)) >= STARTER^MINIMUM^ADJACENT^ALIGN) { 

if(best_new_unmarked_entries_filled < new_unmarked_entriesjilled){ 
best_new_unmarked_entriesj511ed = 

new_unmarked_entries_filled; 

} 

if(DEBUGJSTARTER) fprintf(stdout, "STARTER ALIGN: 
new_unmarked_entries_filled=%d\n , \new_uiimarked_entries_filled); 
epothilonelen = strlen(epothilone.monomersequence); 

for(ii=0; ii< epothilonelen; ii++){ 
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if(DEBUG_STARTER) fprintf(stdout, "STARTER ALIGN:found 

a best alignment 

between epo.monomer[%d]=%c in library[%d].name=%s\n" 3 

ii,epothilone.monomersequence[ii],kk,epothUone.alignedPKSnaine[ii]); 

<• } 

} 

library[kk].recursion_tagged = TRUE; 
^>rintf(stdout,"ALIGN_TARGET: \n"); 
dump_STARTER_aUgn(epothilone^wildcard,wildcards); 
§)rintf(stdout,"ALIGN_TARGET: \n"); 
if(best_new_unmarked entries filled <= 1) { 

§>rintf(stdout/'ALIGN_TARGET: PROBLEM 
bestjnew_unmarked_entxies_filled = %d\n" J best_new_unmarked_entries_filled); 

fprintf(stdout,"ALIGN_TARGET: PROBLEM skipping this STARTER 
entry for library[%d].name=%s\n", 

kk,library[kkj .name); 
library[kk].recursion_tagged = FALSE; 
kk-H-; 
continue; 

} 

/* - fill in the gaps from the library - */ 

/* generate a fresh copy of epothilone in epotemp */ 
epothilonelen ~ strlen(epothilone.monomersequence); 
nfilledmax = strlen(epothilone.monomersequence); 
fprintf(stdout, "nfilledmax=%d\n'\ nfilledmax); 
reset_epotemp(&epotemp,epothiione); 
nfilled = 0; 

for(ii=0; ii< epothilonelen; ii-H-) { if(epotemp.marked[ii] = TRUE) nfilled-H-; } 
if(DEBUG_STARTER) ftrint^stdout/'nfilled from STARTER=%d\n",nfilled); 
for(mmpass = 0; mmpass < nlib; mmpass++) { 

if(mmpass — kk) { continue; } 

reset_epotemp(&epotemp,epothilone); 

nfilled = 0; 

foi^ii^O; ii< epothilonelen; ii++) { 
if(epotemp.marked[ii] = TRUE) nfilled-H-; } 
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if(nfilled >= nfilledmax) { 

output_fresh_alignment(&epotemp); 

} else { 

currentnmarked = nfilled; 

previous__nmarked = nfilled; 

smallest_acceptable_piece = 1; 

while((new_unm arked entries filled = 
maximal_adjacent_alignment(&epotemp,nw 
acceptable_piece)) >= MINIMUM ADJACENT ALIGN) { 

currentnmarked += 

new_unmarked_entries_filled; 

if(DEBUG_MATCH) fprintf(stdout, 
"main: recursion_level=%d, mmpass=%d, previous_nmarked=%d, 
current jnmarked=%d\n M , 

recursion_counter,mmpass,previous_nmarked, currentnmarked); 

} 

nfilled = 0; 

for(ii=0; ii< epothilonelen; ii++) { 
if(epotemp.marked[ii] = TRUE) nfilled++; } 

if(nfilled >= nfilledmax) { 

output_fresh_alignment(&epotemp); 

continue; /* no need to recurse */ 

} 

if(DEBUG_MATCH) fprintf(stdout, "main: 
about to RECURSE: mmpass=%d\n",mmpass); 

library[mmpass].recursion_tagged = TRUE; 

recursion_counterH-; 
reciu-se_through_thejibrary(nfilledmax,epotemp,&epothilone,nwild 
nlib,library,&recursion_counter); 

library[mmpass].reciirsion_tagged = FALSE; 

recursion_counter~; 

> 

library[kk].recursion_tagged = FALSE; 
kk++; 
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}/*nlib*/ 
}/*main*/ 

int recurse_through_the_library( 
int nfilledmax. 
LIB epotemp, 
LIB *epothilone, 
int nwildcard, 

char wildcards[MAXA\TLD][MAXLEN], 

int nlib, 

LIB *library, 

int *recursion_counter) 

{ 

int ii=0; 

int ecount=0,elen=0; 
int mmpass=0; 
int nfilled=0; 
int lcount=0; 

int previous_nmarked=0, cuirent_nmarked==0; 
int smallestacceptablejpiece = 0; 
char *eptr> 
char *clibptr, 

char boundary[MAXNAMELEN]; 
int new_urtmarked_entries_filled=0; 
LIB epotemp_temp; 

ifl[DEBUG_MATCH) fyrintf(stdout, ff RECURSE: recursion_counter=%d, 
n]ib=%d\n",*recursion_connter,nlib); 

elen = strlen(epotemp.monomersequence); 

nfilled = 0; 

for(ii=0; ii< elen; ii++) { if(epotemp.marked[ii] = TRUE) nfilled++; } 

previous_nmarked = nfilled; 

currentnmarked = nfilled; 

smallest_acceptablej>iece = 1; 

if(nfiUed >= nfilledmax) { return 1; } 

for (mmpass = 0; mmpass < nlib; mmpass-H-) { 

if(*recursion_counter>= RECTOR { 
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return 1: 

} 

if(DEBUG_MATCH) ^rint^stdout/RECURSE: 
recursion_counter=%d, mmpass^ 0 d\n T ^*recirrsion_coimter,mmpass); 
if{library[mmpass].recursionJagged = TRUE) { 

if(DEBUG__MATCH) ^rintf^stdout^RECURSE: 
library[%d].recursion_tagged=TRUE; skipping\n",minpass); 
continue; 

reset_epotemp(&epotemp_temp,epotemp); 
elen = strlen(epotemp_temp.monomersequence); 
nfilled = 0; 

for(ii=0; ii< elen; ii++) { if(epotemp_temp.marked[ii] = 
TRUE)nfilled-H-;} 

previousnmarked = nfilled; 
current_nmarked = nfilled; 
while((new ^unmarked _entriesjalled = 
maximal adjacent ali gnment(&epotemp_temp ? nwildcard, wildcards,library, 
mmpass,smallest_acceptable_piece)) >= 1) { 

currcntnmarked +== new_unmarkcd_entries_filled; 
if(DEBUG_MATCH) fprintf(stdout, "RECURSE: recursion Jevel=%d, 
mmpass=%d, previous_nmarked= : %d, current_nmarked=%dVn", 

*recursion_counter J nimpass 9 previous_nmarked, cvurent_nmariced); 
} 

elen = strlen(epotemp_temp.monomersequence); 
nfilled = 0; 

for(ii=0; ii< elen; ii++) { if(epotemp_temp.marked[ii] = 
TRUE)nfilled++;} 

if(nfilled >= nfilledmax) { 

output_fresh_atigrnnent(&epotemp_temp); 

continue; 

} 

library[mmpass].recursion_tagged = TRUE; 
(*recursion_counter)-H-; 
reciu-se_through_meJtibrary(i^^ 
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rds,nlib,library,recursion_counter); 

library[mmpass].recursion_tagged = FALSE; 
(*recursion_counter)~; 
}/* mmpass */ 
} /*recurse_through_the library*/ 
/* 

PURPOSE: 

INPUT: 

OUTPUT: 

returns the size of the largest maximal adjacent set of 
monomers inserted. 

PROCEDURE: 

*/. 

int maximal_adjacent_alignment( 
LIB *epothilone, 
int nwildcard, 

char wildcards[MAXWILJ)][MAXLEN], 
LIB *library, 
int ilib, 

int smallest_acceptable_piece) 
{ 

int ii=0,jj=0,kk=0; 

int ecount=0,lcount=0; 

int epothilonelen=0; 

int nlargestpiece=0,tnlargestpiece=0; 

int hold_this_lcount=0, hold_this_ecount=0; 

int wi!dcardmatch=FALSE; 

char *wptr; 

char *1argestpiece_sptr,*largestpiece_eptr; 
char *hold_thisjplace_eptr, *holdjhis_place_sptr; 
int Iargestpiece_lcount=0 ) largestpiece_ecount=0; 
char *sptr, *eptr, *lptr,*bufptr; 
if(DEBUGJWTLDCARD) { 
if(nwildcard > 0) { 

fprintf(stdout, M maximal_: 
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wildcards[0]=%s\n" 5 wildcards[0]); 
} 

} 

iflPEBUG_ALIGN) fprin 
smallest_acceptablejpiece=%d\n",smallest_acceptable_piece); 
sptr = libraryfilib].monomersequence; 
eptr = epothilone->monomersequence; 
ecount=0; 
lcount=0; 
nlargestpiece=0; 
tnlargestpiece=0; 
hold_this_place_eptr = eptr; 
holdthisecount = ecount; 
while (*eptr!=^ f ){ 

sptr = library[ilib].monomersequence; 

lcount = 0; 

hold_this_place_sptr = sptr; 
hoi d_this_l count = lcount; 
wildcardmatch = FALSE; 
while(*sptr N \0') { 

wildcardmatch = FALSE; 
if(epothilone->marked[ecount] — FALSE) { 

/* code for wildcards added MAS 05-16-00 */ 
wptr = ""; 

if(*eptr = 'X') { wptr = wildcards[0]; } 
else if(*eptr = V) { wptr = wildcards[l]; } 
else if(*eptr = 'Z') { wptr = wildcards[2]; } 
while(*wptr != f \0') { 

if(*wptr = *sptr){ 

wildcardmatch = TRUE; 

break; 

} 

wptrH-; 

} 

if((wildcardmatch = TRUE) || (*eptr = 
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*sptr)) { 

tnlargestpiece++; 

if(DEBUG_ALIGN) fprintf(stdout, 
tr FOUND a match: len=%d, epo(%d, %c), lib[%d].name=%s (%d, %c)\n", 

tnlargestpiece, 

ecount, *eptr, ilib^ibraryplib] Jiame,lcount, *sptr); 

if(tnlargestpiece > nlargestpiece) { 

nlargestpiece = tnlargestpiece; 
largestpiecesptr = 



hold_thisj3lace_sptr; 
holdjhislcount; 
hold_this_place_eptr; 
holdthisecount; 



largestpiecejcount = 
largestpiece_eptr = 
largestpiece ecount = 



if(DEBUG_ALIGN) 
fprintf(stdout, "FOUND a largest piece: len=%d, epo(%d, %c), 
lib[%d].narae=%s (%d, %c)\n", 

nlargestpiece, 

largestpiece_ecount, *largestpiece_eptr, 
ilib,library[ilib].name > largestpiece lcount, *largestpiece_sptr); 

} 

sptr-H-; 
lcount++; 
eptr-H-; 
ecount++; 

} else { 

tnlargestpiece = 0; 
sptH-H; 
lcount-H-; 
/♦NEW*/ 

holdthis _j)lace_sptr = sptr; 
hold_this_lcount = lcount; 
eptr = hold_this_place_eptr; 
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ecount = hold_this_ecount; 

} 

} else { 

tnlargestpiece = 0; 
break; 

} 

} 

tnlargestpiece = 0; 
eptr = hold_this_place_eptr + 1 ; 
ecount = hold Jhis_ecount + 1 ; 
hold_thisjplace_eptr = eptr; 
hold this ecount = ecount; 

} 

if(DEBUG_ALIGN) fprintffctdout/'ALTGN: largest piece match is %d monomers 
from %s\n"^tlargestpiece,library[ihb].name); 

if(DEBUG_ALIGN) fprintf(stdout 5 M ALIGN: largestpiece_ecount=%d, 
largestpiecelcount=%d\n H , 

largestpiece_ecount,largestpiece_lcount); 
if(nlargestpiece >= smallest_acceptable_piece) { 

if(DEBUG_ALIGN) ^rintf(stdout,"ALIGN: incorporated^"); 
lcount = largestpiece_lcount; 
ecount = largestpiece_ecount; 

while(ecount < (nlargestpiece + largestpiece__ecount)) { 
epottalone->alignedsequence[ecount] = 
library[ihb].monomersequence[lcount]; 

strcpy(epothilone->alignedPKSname[ecount],Ubrary[ilib].name); 
strcpy(epothilone->context[eco^^ 

epothilone->marked[ecount] = TRUE; 
if(ecount < (nlargestpiece + largestpieceecount - 
1)) epothilone->boundarytoright[ecount] = FALSE; 
lcount++; 
ecount++; 

} 

} 

return (nlargestpiece); 
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} /*maximal_adj acent_alignment*/ 
/* 

PURPOSE: 

INPUT: 

OUTPUT: 

returns the size of the largest maximal adjacent set of 
monomers inserted. 

PROCEDURE: 

*/ 

int maximal_adjacent_alignment_and_durnp( 
LIB *epothilone, 
int nwildcard, 

char wildcards[MAX WILD] [MAXLENJ, 
LIB library, 
int ilib, 

int smallest_acceptable_piece) 
{ 

int ii=0,jj=0,kk=0; 

int ecount=0,lcount=0; 

int elen=0; 

int epothilonelen=0; 

int nlargestpieceF=0,tolargestpiece=0; 

int hold_this_lcount=0 5 hold_this_ecount=0; 

int wildcardmatch=FALSE; 

ch^r *wptr; 

char *largestpiecesptr,*Iargestpiece_eptr; 

char *holdJhis_place_eptr, *hold this_place_sptr; 

int largestpiece_lcount=0,largestpiece_ecount=0; 

char *sptr, *eptr, *lptr,*bufptr; 

if(DEBUG_WILDCARD) { 
if(nwildcard > 0) { 

fprintfCstdout/'maximalj 
wildcards[0]=%s\n M ,wildcards[0]); 
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} 

} 

fprmtf^stdout/'maximal^ 
smallest_accq)tablej3iece=%d\n ,, J smallestaccq?table_piece); 
sptr = library[ilib].monomersequence; 
eptr = epothilone->monomersequence; 
elen = strlen(epotWlone->monomersequence); 
ecount=0; 
lcount=0; 
nlargestpiece=0; 
tnlargestpiece=0; 
hold_this_place_eptr = eptr; 
holdthisecount = ecount; 
while (*eptr l~Vt) { 

sptr = library[ilib].monomersequerice; 

lcount = 0; 

/* NEW */ 

hold_thisjplace_sptr = sptr; 
holdthislcount = lcount; 
wildcardmatch = FALSE; 
while(*sptr != W) { 

wildcardmatch = FALSE; 
if(epothilone->marked[ecount] = FALSE) { 

/* code for wildcards added MAS 05-16-00 */ 
wptr=""; 

if(*eptr — 'X') { wptr = wildcards[0]; } 
else ifi(*eptr = T') { wptr = wildcards[l]; } 
else if(*eptr = '7!) { wptr = wildcards[2]; } 
while(*wptr!= , \0 f ) { 

if(*wptr = *sptr){ 

wildcardmatch = TRUE; 

break; 

} 

wptrH-; 

} 



WO 01/92991 PCT/US01/17352 

54 

if((wildcardmatch = TRUE) || (*eptr = 

*sptr)) { 

tnlargestpiece-H-; 

ifl(DEBUG_ALIGN) fprintf(stdout, 
"FOUND a match: len=%d, epo(%d, %c), lib[%d].name=%s (%d, %c)\n", 

tnlargestpiece, 

ecount, *eptr, ilib,library[ilib].name,lcouiit, *sptr); 

ififtnlargestpiece > nlargestpiece) { 

nlargestpiece = tnlargestpiece; 
largestpiece__sptr = 

hold_this_placejsptr; 

largestpiece lcount = 

hold this lcount; 

largestpiece__eptr = 

hold_this__place_eptr; 

largestpieceecount = 

hold_this_ecount; 

if(DEBUG_ALTGN) 
fprintf(stdout, "FOUND a largest piece: len=%d, epo(%d, %c), lib(%d, %c)\n", 

nlargestpiece, 

largestpiece_ecount, *largestpiece_eptr, largestpiecelcount, 
*largestpiece_sptr); 

} 

sptrH-; 
lcount-H-; 
eptr++; 
ecount-H-; 

} else { 

tnlargestpiece = 0; 

sptrH-; 

lcount-H-; 

hold_this_j>lace_sptr = sptr; 
holdjhis Jcount = lcount; 
/* NEW */ 

eptr = hold_this_place_eptr; 
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ecount = hold_this_ecount; 

} 

} else { 

tnlargestpiece = 0; 
break; 

} 

} 

tnlargestpiece = 0; 
eptr = hold__this_place_eptr + 1 ; 
ecount = hold_this_ecount + 1 ; 
hold_this_place_eptr = eptr; 
hold_this_ecount = ecount; 
if(DEBUG_ALIGN) { 

fprintf(stdout, M incrementing 
hold_this_place_eptr=%c, 

hold_this_ecount=%d\n",*holdJhisjplace_eptr 9 hold_this_ecount); 
} 

} 

if(DEBUG_ALIGN) §)rintf(stdout,"ALIGN: largest piece match is %d monomers 
from %sW\nlargestpiece,libraiy[ilib].name); 

if(DEBUG_ALIGN) fprint^stdout," ALIGN: largestpiece_ecount=%d, 
largestpiece_lcount=%d\n ,, > 

largestpiece_ecount,largestpiece_lconnt); 
if(nlargestpiece >= smallest__acceptable_piece) { 

if(DEBUG_ALIGN) fprintfCstdout, n ALIGN: incorporated^"); 
lcount = largestpiece_lcount; 
ecount = largestpiece_ecount; 
^)rintf(stdout 5 "ALIGN_TARGET: "); 
for(ii=0; ii<largestpiece_ecount; ii-H-) { 
frrintf(stdout, " "); 

} 

while(ecount < (nlargestpiece + largestpiece_ecount)) { 
epothilone->alignedsequence[ecount] ~ 
librarytilib].monomersequence[lcount]; 

s1rcpy(epotMlone->a]ignedPKSname[^ 
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strcpy(epothilone->context[ecount],library(ilib].context[lcount]); 

epothilone->marked[ecount] = TRUE; 

if(ecount < (nlargestpiece + largestpiece_ecount - 
1)) epothilone->boundarytoright[ecount] = FALSE; 

fyrintf(stdout/Toc , ^library[ilib].monomersequence[lcount]); 

lcount++; 

ecount++; 

} 

for(ii=ecount; ii<elen; ii++) { 
fprintf(stdout," "); 

} 

^printf(stdout,"%s\D n ,library[ilib].name); 

} 

return (nlargestpiece); 
}/*maximal_adjacent_alignrnent_and_diimp*/ 
int output_fresh_alignrnent( 
LIB *epotemp) 
{ 

int ecount=0; 
char *eptr; 

char boundary [MAXNAMELEN]; 

eptr = epotemp->monomersequence; 
ecount - 0; 

epotemp->nboundary = 0; 

strcpy(boundary,epotemp->alignedPKSnarne[ecount]); 
while(*eptr !- \Q<) { 

if(epotemp->boundarytoright[ecount] = TRUE) { 
epotemp->nboundary++; 

} 

ecount++; 
eptr-H-; 

} 

if(epotemp->nboundary > NBOUNDARYCUTOFF) return 1; 
eptr = epotemp->monomersequence; 
ecount = 0; 



WO 01/92991 PCT/US01/17352 

57 

fprintf(stdout s "HIT "); 
while(*eptr != \0>) { 

if(epotemp->alignedPKSname[ecount][0] = W) { 

if(epotemp->boundarytoright[eco\uit] = TRUE) { 
fprintftstdout/'rocrTARGCyos)! 
T \*eptr,epotemp->context[ecount]); 
} else { 

fprintf(stdout 3 "%c:TARG(%s) 
' f ,*eptr,epotemp->context[ecount]); 

} 

} else { 

if(epotemp->boundarytoright[ecount] = TRUE) { 
fprintf(stdout, f, %c:%4s(%s)| 
, ^*eptr,epotemp->alignedPKSname[ecount],epotemp->context[ecount]); 
} else { 

^>rintf(stdout, ,, %c:%4s(%s) 
'^^eptrjepotemp^alignedPKSnametecountJjepotemp^contextfecount]); 

} 

} 

ecount-H-; 
eptr+4-; 

} 

rprintftstdouV'npiece %d",epotemp->nboundary); 

fprint^stdout/Xn"); 

return 1; 
} /*output J5resh_alignment*/ 
int get_library( 
char *libraryfile, 
LIB ^library) 

{ 

int ii=0jj=0Jdk==0,lcount=0; 

int nlib=0; 

char *bui^tr,buf[MAXBUF]; 

char tmonomersequcnce[MAXNAMELEN] ; 

char *lptr,*lptr; 
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FILE *libp; 

for(kk==0; kk < MAXLIBJENTRIES; kk-H-) { 
library[kk].recursion_tagged = FALSE; 
for(ii=0; ii<MAXNAMELEN; ii-H-) { 
library[kk].name[ii] = W; 
library[kk].monomersequence[ii] = *\0 f ; 
library[kk].alignedsequencefii] = W; 
for(jj=0; jj<MAXLEN; { 

Iibrary[kk].aIigaedPKSnameOj][ii] = W; 
library[kk].niarked[jj] - FALSE; 
library[kk].context[jj][0] = W; 
library[kk].contextQj][l] = *\0'; 
Hbrary[kk]xontext[ij][2] = W; 
library[kk].contextQj][3] = W; 

} 

} 

library[kk] .nboundary = 0; 

/* read in the library from PKS.lib */ 

if(NULL=^ibp===fopen(libraryfile, , 'r ,, ))) { 

fprintf(stdout,''TRY AGAIN; couldn't open %s\n",libraryfile); 

nlib=0; 

exitO; 

} 

nlib=0; 

while(nUb < MAXLBENTRIES) { 

ifi^nJLL=fgets(buf } skeof(biif),libp)) break; 
burbtr = buf; 

if(*bufjptr = "#') continue; 
if(*bufptr!=V){ 

lptr = library[nlib].name; 

wbile((*bufptr !=' ')&& (*bufptr != W) && 

(*bufptr 1= W)){ 

*lptr++ = *bufptr++; 

} 
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*lptr = '\0*; 

if((*bufptr != \W) && (*bufptr != V)) bufptr++; 
lptr = library[nlib].monomersequence; 
while((*bufptr N ' ") && (*buiptr != \0 r ) && 

(*bufbtr != V)){ 

/* This code specifically deletes inter-modular 
double bonds, optional.*/ 

if(*bufptr !='=') { 

*lptr++ = *bufptr-H-; 

} else { 

bufptr++; 

} 

} 

*lptr='\0'; 

ilX(*bufptr != \V) && (*bu$tr != ^n 1 )) bufptrH-; 
lptr = library[nlib].annotatedsequence; 
while((*bufrtr != ' •) && (*bufptr != "NO') && 

(*bufptr !=V)){ 

*lptr++ = *bufptr++; 

} 

*lptr = W; . 

if((*bufptr != W) && (*bufptr != V)) bufptrH-; 

fprintf(stdout,"LIBRARY(%d) %s: 
%s\n"^b,library[nlib].name,library[nlib].monomersequence); 

fprmtf(stdouCLBRARY(%d) %s: 
%s\n"^uib s library[nlib].nanie,library[nlib].annotatedsequence); 

nlib++; 

} 

} 

fcloseflibp); 

for(kk=0; kk< nlib; kk++) { 
Icount = 0; 

lptr = library[kk].monomersequence; 
while(*lptr!='\0 , ){ 

if(lcount = 0) { 
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Iibrary[kk].coatext[lcount][0] = 
library[kk].cont6Xt[lcount][l] = *lptr; 
library[kk]xontext[lcount][2] = *(lptr + 1); 

} else { 

if(lcount = 
(strlen(library[kk] .monomersequence) - 1 )) { 

library [kk].context[lcount][0] = 



*(lptr-l); 
*lptr; 



*(lptr-l); 
*lptr, 

*(lptr+l); 



Iibrary[kk].context[lcount][l] = 
library[kk].context[lcount][2] = 

} else { 

library [kk]. context [lcount][0] = 
librarytkk].context[lcount][l] = 
library[kk].context[lcount][2] = 

} 



} 

library[kk].context[lcount][3] = f \0'; 

lptrf+; 

lcountf-+; 

} 

fprintffstdout/'LBRARY^/od) %s: 
%sW'Jkk,Hbrai7[kk].name > libraTy[kk].monomerseque!nce); 

fprintf(stdo<TmRARY(%d) %s: 
%sW'4dfc,library[kk].name,library[kk].anriotatedsequence); 

for(jj=0; jj< strlen(librar3^[kk].monomersequence); { 
^rintf(stdout;'(%s) M ,Iibrary[kk].context[jj]); 

} 

^rintfltstdout.-Nn"); 

} 

return nlib; 
}/*get_library*/ 
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int dump_STARTER_align( 
LIB epothilone, 
int nwildcard, 

char wiIdcards[MAXWILD]pvIAXLEN] 

) 

{ 

int elen=0; 

int ecount^hold^ecount^O; 
int wildcardmatch=FALSE; 
char *sptr,*eptr,*wptr; 

elen = strlen(epothilone.monomersequence); 
eptr = epothilone .monomersequence; 
fprintfCstdout, "ALIGN_TARGET: 
while(*eptr != *\0') { 

fyrintf(stdout; , %c n ,*eptr); 
eptr++; 

} 

fyrintf(stdout,"\n"); 
ecount=0; 

eptr = epothilone.monomersequence; 
sptr = epothilone.alignedsequence; 
Q)rintf(stdout, ,r ALIGN_TARGET: 
while(*eptr != '\0') { 

wildcardmatch = FALSE; 

wptr = nn ; 

if(*eptr = 'X') { wptr = wildcards[0]; } 
else if(*eptr = 'Y') { wptr = wildcardsfl]; } 
else if(*eptr = 'Z') { wptr = wildcards]?]; } 
while(*wptr { 

ifl[*wptr = *sptr){ \ 

wildcardmatch = TRUE; 

break; 

} 

wptr++; 

} 
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if(wildcardmatch = TRUE) { 
Q)rintf(stdout, ,, : ,r ); 

} else { 

if(*eptr = *sptr) { 

Q}rintf(stdout, ,, |"); 

} else { 

frrintf(stdout, M "); 

} 

} 

eptrH-; 
sptrH-; 

} 

fprintf(stdout, rt \n M ); 

fprintf(stdout, "ALIGNJTARGET: "); 
eptr = epothilone.monomersequence; 
ecount=0; 

while(*eptr != \W) { 

if(epothilone.alignedsequence[ecount] = W) { 
fyrintf(stdout, n "); 

} else { 

fprintf(stdout/ , %c n > epothilone.aIignedsequence[ecount]); 
hold ecount = ecount; 

} 

eptrf+; 
ecount-H-; 

} 

$rintf(stdout," 

%s\n",epothilone.alignedPKSnamethold_ecount]); 

Q)rintf(stdout; , STARTER_ALIGN:\n M ); 
}/Mump_STARTER_align*/ 

reset_epotemp( 
LIB *epotemp, 
LIB epothilone) 

{ 

int jj=0, elen=0; 
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elen = strlen(epothilone.monomersequence); 
strcpy(epotemp->name,epothiloiie.name); 

strcpy(epotemp->monomersequence 5 epothilone.monomersequence); 
strq)y(epotemp->alignedsequence,epotiiilone.alignedsequence); 
epotemp->nboundary = 0; 
for(jj=0; jj< elen; jj-H-) { 
strcpy(epotemp->ahgnedPKSname[jj],epotta^^ 

epotemp->marked[jj] = epothilone.marked[jj]; 

strq)y(epotemp->context|jj] 5 epotMlone.context[ij]); 

qpotemp->boundarytoright(jj] = epothilone.boundarytoright[jj]; 

} 

} /*reset_epotemp*/ 

EXAMPLE 6 

Source Code: 
#include <stdio.h> 

/* ~siani/programs/morph/morph4.c 

PURPOSE: To recursively traverse all the entries in PKS.lib, generating all feasible 
combinations of PKS modules to make the TARGET (e.g., epothilone). 

INPUT: -b number_boundary_cutoff: lets user set the maximum number of 
boundaries in output lines. This defaults to 5 (#define NBOUNDARY CUTOFF 5) which is a 
reasonable assumption for something of the length of epothilone (8 modules). However, when 
looking at disco-dennolide which has 1 1 modules, a cutoff of 5 sometimes results in too few 
output lines; it is too restrictive. 

-d allows one to ignore the inter-modular doublebonds in the library file. 

-1 libraryfile: tab-delimited CHUCKLES-coded polyketides file with the following 
columns 

1 . polyketide name 

2. plain CHUCKLES 

3. annotated CHUCKLES (contains information about post-synthetic 

modifications) 

4. source organism 

-n targetname: user-defined name (e.g., epoD) 

-t targetsequence: CHUCKLES-coded polyketide of desired TARGET (e.g., 
MEMLJDGE) 

-w, -x, -y, -z sets of wildcards: sets of monomers for particular positions appearing in 
targetsequence. The wildcards can effectively be used for analoging the TARGET polyketide. 

Hard-coded parameters which may be reset (requires recompiling): 
#defineNBOUNDARY_CUTOFF 5 

NB OUND AR Y_CUTOFF determines the maximum number of non-native inter- 
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modular interfaces which are contained in the output. This is now set to 5, but may be increased 
when the user does not care about inefficiencies introduced by these interfaces or when the 
targetsequence is very lengthy. 

#defineRE(^SION_COUNTER_CUTOFF 2 

RECURSION_COUNTER_CUTOFF specifies the number of levels of recursion 
(defaults to 0, 1, 2) acceptable for the run. This limit must be set since the large PKS library can 
result in recursion that will combinatoriaUy explode. Because of the multi-directionality of the 
alignments (using every library entry as a STARTER), there is no need to go beyond 2 levels of 
recursion. However, there may be cases in the future where this number should be increased. 
Note that while recursion will eventually terminate without this parameter, runs with a library 
over about 20 PKS entries may run for years on a reasonably fast computer. 

OUTPUT: All combinations of modules that meet parameters set by user. 

Example output from MEMLJDGE (epothilone D) using subset of PKS.lib. 
Vertical bars indicate non-native inter-modular interfaces. Last column contains the number of 
"pieces" that are needed to put together the PKS. 

Names of PKSs have been abbreviated to fit them in these comments. 
HIT M:3atyl(FMN)| E:tedan(GEH)| M:aldga(BML) L:aldga(MLG)| J:aldga(GJD) 
D:aldga(JDL)| G:tedan(JGE) E:tedan(GEH)| 5 

HIT MralbMl(LME) E:albMl(MEJ)| M:albMl(LME)| L:aIdga(MLG)| J:aldga(GJD) 
D:aldga(JDL)( G:tedan(JGE) E:tedan(GEH)| 5 

HIT M:albMl(LME) E:albMI(MEJ)| M:aIdga(BML) L:aIdga(MLG)| . J:aldga(GJD) 
D:aldga(JDL)| G:3atyl(NGO)| E:albMl(MEJ)| 5 

HIT M:albMl(LME) E:albMl(MEJ)| M:aldga(BML) L:aldga(MLG)| J:aldga(GJD) 
D:aldga(JDL)| G:aldga(LGJ)| E:albMl(MEJ)| 5 

HIT M:albMl(LME) E:aJbMl(MEJ)| M:aldga(BML) L:aldga(MLG)| J:aldga(GJD) 
D:aldga(JDL)| G:aldga(LGJ)| E:albMl(MEJ)| 5 

USAGE: 

morph3 -1 libraryfile -n targetname -t targetsequence [-w W-wildcards] [-x X- 
wildcards] [-y Y-wildcards] [-z Z-wildcaxds] -d 

examples: 

# generate combinations that yield epothilone D 

%morph3 -1 PKS.lib -n epoD -t MEMLJDGE > omorph3_epoD 

%egrep HIT omorph3_epoD | sort | uniq | sort +10 -1 1 > 
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omorph3_epoD.uniq.sort 

%egrep ALIGN JTARGET omorph3jspoD > 
omorph3_epoD_STARTER_ALIGN 

# generate combinations that yield epothilone D with a CI 3-hydroxyl 
%moiph3 -1 PKS.lib -n epoD-130H -t MEXLJDGE -x ABCD > oepoD-130H 
%egrep HIT oepoD-130H | sort | uniq | sort +10 -1 1 > oepoD-130H.uniq.sort 
%egrep ALIGN_TARGET oepoD-130H > oepoD- 1 30H_STARTER_ALIGN 

# generate combination that yield epothilone with the following wildcards (set 1) 
%morph3 -1 PKS.lib -n epoD-setl -t MEXYZDgE -x ABCD -y LEFESF -z 

JACGM > oepoD-setl 

%grep HIT oepoD-setl | sort | uniq | sort +10 -1 1 > oepoD-setl.uniq.sort 

# generate combination that yield epothilone with the following wildcards (set 2) 
%morph3 -1 PKS.lib -n epoD-set2 -t MEXYZDgE -x JK -y EF -z JACGM > 

oepoD-set2 

%grep HIT oepoD-set2 | sort | uniq | sort +10 -1 1 > oepoD-set2.uniq.sort 
LIMITATIONS: 

Current implementation cannot handle intra-modular modifications/splitting 
because morph is operating at the monomer level. Future implementations could convert the 
CHUCKLES-encoded strings into the corresponding and equivalent SMILES and then perform 
more complex chemical analysis of the PKS molecular graphs. Currently, inter-modular double 
bonds are present in the library, but are ignored by the morph program. 

MODIFICATIONS: 

+ added ability to include user-defined wildcards (X, Y, or Z) on the 

command line. MAS 05-1 6-00. 
+ added additional wildcard (W). MAS 05-30-00. 
+ added addition (summary) column to HIT output list. MAS 05-30-00. 
+ added command line argument for suppressing the inter-modular double bonds 
in the library. Default is not to treat these as separate modules. MAS 05-3 1-00. 

+ added column that contains the length of the largest matching fragment MAS 

06-05-00 
*/ 
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#defineTRUE 1 
#defineFALSE 0 
#defineDEBUG_MATCH FALSE 
#defineDEBUG_STARTER FALSE 
#defineDEBUG_ALIGN FALSE 
#defineDEBUG_RECURSE FALSE 
#defineDEBUG_WILDCARD FALSE 

#defineMAXLEN 80 
#defineMAX_TYL_LEN 6 
#defineMAX_EPO_LEN 6 
#defineMAXNAMELEN 500 
#defineMAX_LIB_ENTRIES 500 
#defmeMAXWILD 4 
#defineMAXBUF 1000 



#defineNBOUNDARY_CUTOFF 5 
#defineRECURSION_COUNTER_CUTOFF 2 
#defmeSTARTER_MIMMUM_ADJACENT_ALIGN 2 
#defmeMINIMUM_ADJACEMT_ALIGN 2 



typedef struct lib { 

char name[M AXNAMELEN] ; 

char monomers equence[MAXNAMELEN] ; 

char annotatedsequence[MAXNAMELEN]; 

char alignedsequence[MAXNAMELEN]; 

char ahgnedPKSname[MAXLEN][MAXNAMELEN]; 

int boimdarytorightlMAXNAMELEN]; 

int markedfMAXLEN]; 

char contextpviAXLEN][4]; 

int recursionjagged; 

int nboundary; 

}LEB; 
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main(int argc, char **argv) 
{ 

int ii=0, jj=0, kk=0, 11=0; 

int nlib=0; 

int ecount=0; 

int nfi lled=0,nfi lledmax=0; 

int epothilonelen=0; 

int nlargestpiece=O 3 talargestpiece==0; 

int mmpass=0; 

int lcount=0; 

int new_unmarked_entries_fiiled = =0; 

int recursioncounter = 0; 

int nwildcard=0; 

int best_new_unmarked_entries_filled = 0; 

int smallest_acceptablej)iece - 0; 

int current_nmarked=0,previous_nmarked=0; 

int inter_modular_db_flag_off = FALSE; 

int nboundary_cutoff^NBOUNDARY_CUTOFF; 

char *sptr, *eptr, *lptr,*bufptr; 

char *clibptr; 

char *libraryfile; 

char *targetsequence,*targetname; 

char buflMAXBTJF]; 

char wildcards[MAXWILD][MAXLEN] ; 

FILE *libp; 

LIB epotemp; 

LIB Ubrary[MAX_LIB__ENTRIES]; 

LIB epothilone; 



char *progname; 

char **filelist, **fileptr, 
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libraryfiie = ""; 
targetsequence = ""; 
targetname = ""; 

for(ii=0; ii<MAXWILD; ii++) { 

foT(y=0;jj<MAXLEN;jj++) { 
wildcards[ii][}j] = \0'; 

} 

} 

/* process arguments */ 

filelist = fileptr = (char **)(malloc(argc * sizeof(*argv))); 

progname = *argv++; . 
if{argc<2) { 

fprintf(stderr,"usage:%s [-b nboundary_cutoff] [~d] -1 libraryfile -n targetname -t 
targetsequence [-w W-wildcards] [-x X-wildcards] [-y Y-wildcards] [-z Z-wildcards] 
\n",progname); 

exitO; 

} 

while(argc- > 1) { 

if(argv[0][0] = •-• && argv[0][l] 1= \V) { 
/* handle option */ 

*-H-(*argv); /* advance past the minus */ 
switch(**argv) { 

case V: /* get number of boundaries cutoff for output of 

alignments */ . 

argv++; argc-; 

sscanf(argv[0], ,T %d n ,&nboundary_cutofO; 
fprmtf(stderr,"-b: 
nboimdary_cutoff=%d\n" ^boundary cutofi); 

break; 

case f d': /* ignore inter-modular double bonds in the library file */ 
biter_modular_dbJlag_off = TRUE; 
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fprintf(stderr/ -d: inter-modular double bonds ignoredAn"); 
break; 

case T: /* get library input filename (PKS.lib) */ 
argv++; argc--; 
libraryfile = argv[0]; 

fprintf(stderr 5 "-l: libraryfile=%s\n" > libraryfile); 
break; 

case f ri: /* get target name string */ 
argv-H-; argc--; 
targetname = argv[0]; 

fprintf(stderr,"-t: targetname=%s\n M ,targetname); 
break; 

case V: /* get target sequence string */ 
argv-H-; argc-; 
targetsequence = argv[0]; 

fprintffetderr/'-t: targetsequence=%s\n'\targetsequence); 
break; 

case W: /* get a wildcard string */ 
argv++; argc--; 
strcpy(wildcards[0] 5 argv[0]); 

fprintf(stderr,"-w: wildcards[%d]=%s\n" 5 0,wildcardstO]); 

nwildcard++; 

break; 

case *x': /* get a wildcard string */ 
argv-H-; argc-; 
strcpy(wildcards[l],argv[OJ); 

fprintf(stderr,"-x: wildcards[%d]=%s\n",l,wildcards[l]); 

nwildcard++; 

break; 

case y : /* get a wildcard string */ 
argv++; argc-; 
strcpy(wildcards[2],argv[0]); 

fyrintf(stderr/'-y:wildcards[%d]=%s\n",2 5 wildcards[2]); 

nwildcard-H-; 

break; 
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case 'z 1 : /* get a wildcard string */ 
argv-H-; argc~; 
strcpy(wildcards[3],argv[0]); 

fprintf(stderr, M -z: wildcards[%d]=%s\n T, ,3,wildcards[3]); 

nwildcard-H-; 

break; 

case 'W f ; /* get a wildcard string */ 
argv++; argc-; 
strcpy(wildcards[0],argv[0]); 

fyrintf(stderr/'-w: wildcards[%d]=%s\n",0,wildcards[0]); 

nwildcard++; 

break; 

case 'X': /* get a wildcard string */ 
argv++; argc-; 
strcpy(wildcards[l],argv[0]); 

fprintf(stderr,"-x: wildcards[%d]-%s\n n , 1 ,wildcards[ 1 ]); 

nwildcard++; 

break; 

case V: /* get a wildcard string */ 
argv++; argc-; 
strcpy(wildcards[2] ? argv[0]); 
, fprintftstderr/'-y: wildcards[%d]=%s\n",2,wildcards[2]); 
nwildcard++; 
break; 

case 'Z: /* get a wildcard string */ 
argv++; argc-; 
strcpy(wildcards[3],argv[0]); 

fprintf(stderr/'-z: wildcards[%d]=%s\n ,, ,3 > wildcards[3]); 

nwildcard-H-; 

break; 

default: 

jprintf(stdeiT, M %s unknown option; ignored\n",*argv); 

}/*switch*/ 
} else { /* a regular filename */ 
*fileptr-H- = *argv; 
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*fileptr = NULL; 

} 

argv++; 

}/*while*/ 



ifi[nwildcard > 0) { 

forCiiH); ii< nwildcard; ii++) { 

^rintf(stdeir/ , wldcardsf%d]=%s\a",ii,wildcards[ii]); 
j^rintf(stdout;Vildcards[%d]=%s\h M 5 ii,wildcards[ii]); 

} 

} 

epothilone.nboundary = 0; 

for(ii=0; ii<MAXNAMELEN; ii-H-) { 

epothilone.name[ii] = r \0 ? ; 

epothilone.monomersequence[ii] = W; 

epothiIone.aIignedsequence[ii] = 

epothilone.boundarytori^it[ii] = TRUE; 

for(jj=0; jj<MAXLEN; { 

q>othilone.alignedPKSname[ij][ii] = \0*; 
epothilone.marked[jj] = FALSE; 
epothilone.context[j}][0] = 'NO'; 
epothilonexontextftj][l] = *\0'; 
epothilonexontext[jj][2] = l \0 ! ; 

} 

} 

strcpy(q>othilone.name,targetname); 
strq)y(epothi]one.monomersequence,targetsequence); 

fprintf(stdout, "TARGET: %s\n", epothilone.monomersequence); 



ecount = 0; 
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eptr = epothilone.monomersequence; 
while(*eptr 

if(ecount = 0) { 

epothilonexontext[ecount][0] = f -'; 

epothilonexontext[ecount][l] = *eptr; 

epothilonexontext[ecount][2] = *(eptr + 1); 

} else { 

if(ecount = (strlen(epothilone.monomersequence) - !)){ 
epothilone.context[ecouiit][0] = *(eptr- 1); 
epothilone.context[ecount] [ 1 ] = *eptr; 
epothilone.context[ecount][2] = 

} else { 

epothilone.context[ecount][0] = *(eptr - 1); 
epothilone.context[ecount][l] = *eptr; 
epothilone.context[ecount][2] = *(eptr-f- 1); 

} 

} 

epothilone.context[ecount][3] = 

eptrH-; 

ecount-H-; 

} 

for(ii=0; ii<ecount; ii-H-) { 

^rint^stdout/X^oS^^epothilone.contextfii]); 

} 



/* library */ 

nlib = get_Ubraiy(Ubraryfile,bT>rary^ 
^rintf^stdout^nlib^/odW^nlib); 



kk=0;while(kk<nlib) { 
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/* Z ero out the epothilone entry with respect to a new alignment */ 

for(ii=0; ii<MAXNAMELEN; ii++) { 

epothilone.alignedsequence[ii] = '\0 r ; 
epothilone.boundarytoright[ii] = TRUE; 
for(ij=0: jj<MAXLEN; jij++) { 

epothilone.alignedPKSname[jj][ii] = W; 
epothilone.marked[jj] = FALSE; 

} 

} 



/* reset the context back to that in epothilone */ 

ecount = 0; 

eptr = epothilonejnonomersequence; 
while(*eptr != '\0'){ 

if(ecount = 0) { ■ 

epothilone.context[ecount][0] = 

epothilone.context[ecount][l] = *eptr; 

epothilone.context[ecount][2] = *(eptr + 1); 

} else { 

if(ecount = (strlen(epothilone.monoraersequence) - !)){ 
epothilone.context[ecount][0] = *(eptr- 1); 
epothilone.context[ecount][l] = *eptr; 
epothilone.context[ecount][2] = ; 

} else { 

epothilone.context[ecount][0] = *(eptr - 1); 
, epothilone.context[ecount][l] = *eptr; 
epothilonexontext[ecount][2] = *(eptr + 1); 

} 

epothilone.context[ecount][3] = 'W; 

eptrf+; 

ecount++; 

} 



/* 



align STARTER (current library entry) and epothilone 
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sptr = library[kk].monomersequence; 
lcount=0; 

wiule(*sptr!=^{ 

^rint^(stdout/library[%d].monome^sequence[%d]=%c\n , ^ 
kk,]count 5 library[kk].nionomersequence[lcount]); 

sptrH-; 
lcount-H-; 

} 

/* Call maximal^adj acent alignment until it no longer returns more than 
two adjacent modules. There is really no reason to try to extract 
individual modules because this will be done as part of the 
recursive filling of spaces from the library. 

*/ . 
smallestacceptable __piece = 2; 
eptr = epothilone.monomersequence; 
fprintf(stdout, "ALIGN_TARGET: "); 
while(*eptr != W) { 

fyrintf(stdout,"%c ,, ,*eptr); 

eptrH-; 

} 

iprintf(stdout,'V'); 

§>rintf(stderr,"aligning %d %s\n M ,kk, library[kk].name); 

best_new_unmarked_entries_filled = 0; 
while((new_iminarked_entries filled = 
maximal_adjacent_alignment_and_dum^ 

acceptable_piece)) >= STARTER_MINIMIM_ADJACENT_ALIGN) { 

if(best_new_unmarked_entries filled < new_unmarked_entriesj511ed){ 
best_new_unmarked_entries_filled = 
newunmarkedentriesfilled; 

if(DEBUG_STARTER) fprintf(stdout, "STARTER ALIGN: 
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new_iramarked_entriesJ51Ie^ 

epothilonelen = strlen(epothilone.monomersequence); 
for(ii=0; ii< epothilonelen; ii++){ 

if(DEBUG_STARTER) fprintf(stdout, "STARTER ALIGN:found 
a best alignment between epo.monomer[%d]=%c in Iibrary[%d] Jiame=%s\n", 

ii,epothilone,monomersequence[ii],ldc,epothilone.alignedPKSname[ii]); 
} 

} 

library[kk].recursion_tagged = TRUE; 

ft)rintf(stdout,"ALIGN_TARGET: \n"); 
dump_STARTER_align(epothilone^iwildcard,wildcards); 
fprintf(stdout } "ALIGN_TARGET: \n"); 



if(best_new_unmarked_entries Jilled <= 1) { 

§)rintf(stdout/ALIGN_TARGET: PROBLEM 
best new_unmarked_entries filled = 0 /od\n",best_new_unniarked_entries_fi]led); 

§>rintf(stdout, r, ALIGN_TARGET: PROBLEM skipping this STARTER 
entry for library[%d].name=%s\n , \ 

ldc,library[kk] .name); 
library[kk] .recursiontagged = FALSE; 
kk++; 
continue; 

} 



/* fill in the gaps from the library */ 

/* generate a fresh copy of epothilone in epotemp */ 
epothilonelen = strlen(epothilone.monomersequence); 
nfilledmax = strlen(epothilone.monomersequence); 
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tyrintftstdout, ,, nfilledmax=%d\n ,, 5 nfilledmax); 

reset_epotemp(&epotemp,epothilone); 
nfilled = 0; 

for(ii=0; ii< epothilonelen; ii-H-) { if{epotemp.marked[ii] = TRUE) nfilled++; } 
ifipEBUGjSTARTER) frrintf(stdout/'nfilled from STARTER=%d\n",nfilled); 

for(mmpass = 0; mmpass < nlib; mmpass-H-) { 

if(mmpass — kk) { continue; } 
resetepotemp(&epotemp,epothilone); 



nfilled = 0; 

for(ii=0; ii< epothilonelen; ii-H-) { if(epotemp.marked[ii] = TRUE) 



nfiUed-H-; } 



infilled >= nfilledmax) { 

output^fresh^alignmentC&epotempjnboundarycutoff); 

} else { 

currentjnmarked = nfilled; 
previous_nmarked = nfilled; 
smallest_acceptablej>iece = 2; 
while((new_unmarked_entries_filled = 
maximal_adjacent_aUgnment(&epoto 
lejiece)) >= IvnNIMUM ADJACENT ALIGN) { 

currentjnmarked += new_unmarked_entries_filled; 
if(DEBUG_MATCH) fprintf(stdout, "main: 
recursion_level=%d, mmpass=%d, previouSjiimarked^d, current_nmarked=%d\n", 

recursion_coimter,mmpass,previous_nmarked, 

currentnmarked); 

} 



WO 01/92991 



77 



PCT/US01/17352 



TRUE) nfilled-H-; } 



nunpass= : %d\n t, 5 mnapass); 



nfilled-0; 

for(ii=0; ii< epothilonelen; ii++) { if(epotemp.marked[ii] = 

if(nfilled >=nfilledmax) { 

outputfresh _aligmnent(&epotemp,nboundary_cutoff); 
continue; /* no need to recurse */ 

} 

if(DEBUG_MATCH) fprintf(stdout, "main: about to RECURSE: 

library[mmpass].recursion_tagged — TRUE; 
recursion_counter++; 



recurse throughjhe librar^^ 
ibrary,&recursion_counter^iboundary cutoff); 

library[mmpass] ,recursion_tagged = FALSE; 
recursion_counter~; 

} 

library [kk].recursion_tagged = FALSE; 
kk++; 
}/*nlib*/ 

}/*main*/ 



PURPOSE: 
INPUT: 
OUTPUT: 
PROCEDURE: 



*/ 
int 



recurse_thix>ugh_the_library( 
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int nfilledmax, 
LIB epotemp, 
LIB *epothilone, 
int nwildcard, 

char wildcards[MAXWILD] [MAXLEN], 

int nlib, 

LIB *library, 

int *recursion_counter, 

int nboundary_cutofi) 

{ 

int ii-0; 

int" ecountN^elen^O; 

int mmpass=0; 

int nfilled=0; 

int lcount=0; 

int previous_nmarked=0, current _nmarked=0; 

int smallest_acceptablejpiece = 0; 

char *eptr; 

char *clibptr; 

char boimdary[MAXNAMELEN]; 

int new_unmarked_entries_filled=0; 

LIB epotemp_temp; 

if(DEBUG_MATCH) j^rintf(stdout, M RECURSE: recursion_counter=%d, 
nlib=%d\n M ,*recursion_counter^ilib); 

elen = strlen(epotempjnonomersequence); 
nfilled = 0; 

for(ii=0; ii< elen; ii++) { if(epotemp.markedf ii] = TRUE) nfilled++; } 
previousjamarked = nfilled; 
currentnmarked = nfilled; 
smallest_acceptable_piece = 1; 
if(nfilled >= nfilledmax) { return 1; } 



for (mmpass = 0; mmpass < nlib; mmpass+-f) { 



WO 01/92991 



79 



PCT/US01/17352 



if(*recursion_coimter>=RECURSION_COUNTER_CUTOFF) { 
return 1; 

} 

if(DEBUG_MATCH) fprintf(stdouV'RECURSE: recursion j;ountep^%d, 
nmipass=%d\n , ^*recursion_coimter,ininpass); 

if(library[mmpass].recursion_tagged = TRUE) { 

if(DEBUG_MATCH) ^rintf(stdout,"RECURSE: 
library[%d] xecursion_tagged=TRUE; skipping\n",mmpass); 
continue; 

} 

reset _epotemp(&epotemp_temp,epotemp); 



elen = strlen(epotemp_temp.monomersequence); 
nfilled = 0; 

for(ii=0; ii< elen; ii++) { if(epotemp Jemp.marked[ii] = TRUE) nfilled++; } 
previous nmarked = nfilled; 
current_nmarked = nfilled; 

while((new_unmarked_entries_filled = 
maximal_adjacent_alignment(&epotemp_temp, nwildcard,wildcards,library, 
nunpass,smaUest_acceptable_piece)) >= 1) { 

current_nmarked += new_unmarked_entries__filled; 
ifiQDEBUG_MATCH) fprintf(stdout, "RECURSE: recursionJevel=%d, 
mmpass=%d, previous_nmarked=%d, current_nmarked=%d\n n , 

, *recursion_counter,mmpass,previous_nmarked, 
current_nmarked); 

} 



elen = strlen(epotemp_temp.monomersequence); 
nfilled = 0; 

for(ii=0; ii< elen; ii++) { if(epotempJemp.marked[ii] = TRUE) nfilled-H-; } 
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if(nfilled >= nfilledmax) { 

outputJfreshaligimient(&epotemp_te^ 
continue; 

} 

library[mmpass] .recursionjagged = TRUE; 
(*recursion_counter)-H-; 

rccurseJhrough_theJibrary(n^^ 
nlib,library,recursion_counter,nboundary_cutoff); 

library[mmpass].recursion_tagged = FALSE; 
(*recursion_counter)~ ; 

}/* mmpass */ 

} /*recurse_through_the_library*/ 

/* 

PURPOSE: 

INPUT: 

OUTPUT: 

returns the size of the largest maximal adjacent set of monomers inserted. 
PROCEDURE: 

*/ 

int maximal_adjacent_alignment( 
LIB *epothilone, 
int nwildcard, 

char wildcardsfMAXWILD] [MAXLEN], 
LIB *library, 
int ilib, 

int smallest_acceptable_piece) 
{ 

int ii=0 Jj=0, kk=0; 
int ecount==0,lcount=0; 
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int epothilonelen=0; 

int nlargestpiece=O J tnlarges^)iece=0; 

int hold_this_lcount==0, hold_this_ecount=0; 

int wildcardmatch=FALSE; 

char *wptr; 

char *largesfpiece_sptr,*largestpiece_eptr; 

char *hold_this_place_eptr, *hold_this_place_sptr; 

int largestpicce^lcount^OJargestpieceecount^; 

char *sptr, *eptr, *lptr,*bufptr; 

if(DEBUG__WILDCARD) { 
if(nwildcard > 0) { 

fprintf(stdouV f maximal_: wildcards[0]=%s\n",wildcards[0]); 

} 

} 

if(DEBUG_ALIGN) fyrintf(stdout,"maximal_adjacent_ahgnment: 
sma11est_acceptable_piece=%d\n w ,smallest acceptable_piece); 

sptr = library[ilib].monomcrsequence; 

eptr = epothilone->monomersequence; 

ecount=0; 

IcounM); 

nlargestpiece=0; 

tnlargestpieceO; 

hold_this__place_eptr = eptr; 

hold_this_ecount = ecount; 

while (*eptr != W) { 

sptr = library[ilib].monomersequence; 
lcount = 0; 

hold__this_place_sptr - sptr; 
hold_thisJcount = lcount; 
wildcardmatch = FALSE; 
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while(*sptr != W) { 

wildcardmatch = FALSE; 
if(epothilone->marked[ecount] — FALSE) { 

/* code for wildcards added MAS 05-16-00 */ 

wptr= m '; 

if(*eptr = 'W) { wptr = wildcards[0]; } 
else if(*eptr = 'X') { wptr = wildcards[l]; } 
else if(*eptr = TT) { wptr = wildcards[2]; } 
else if(*eptr = 'Z') { wptr = wildcards[3]; } 

while(*wptr != '\0*) { 

if(* wptr = *sptr) { 

wildcardmatch « TRUE; 
break; 

} 

wptr++; 

} 

if((wildcardmatch — TRUE) || (*eptr = *sptr)) { 



tnlargestpiece-H-; 

if(DEBUG_ALIGN) fprintf(stdout, "FOUND a match: 
len=%d, epo(%d, %c), lib[%d].name=%s (%d, %c)\n", 

tnlargestpiece, ecount, *eptr, 

ittb,library[ilib].name,lcount, *sptr); 

if(tnlargestpiece > nlargestpiece) { 

nlargestpiece = tnlargestpiece; 

largestpiece_sptr = holdthis__place_sptr; 

largestpiece_lco\rat — hold_this_lcount; 

largestpiece_eptr = hold_tbis_place eptr; 

largestpiece_ecount = holdjhisecount; 

if(DEBUG_ALIGN) fprintf(stdout, "FOUND a 
largest piece: len=%d, epo(%d, %c), lib[%d] Jiame=%s (%d, %c)\n", 

nlargestpiece, largestpieceecount, 
*largestpiece_eptr, ilib,b*brar>'[ilib].name,largestpiece_lcount, *largestpiece_sptr); 

} 
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sptrH-; 
lcount++; 
eptr++; 
, ecount++; 
} else { 

tnlargestpiece = 0; 
sptr++; 
lcount-H-; 
/* NEW */ 

hold_this_pIace_sptr = sptr; 
hold_this_lcount = lcount; 
eptr = hold_this_place_eptr; 
ecount = hold_this_ecount; 

1 } 

} else { 

tnlargestpiece = 0; 
break; 

} 

} 

tnlargestpiece = 0; 
eptr = hold_this_place_eptr + 1; 
ecount = hold_Jhis_ecount + 1 ; 
holdjhis _place_eptr = eptr, 
hold_this_ecount = ecount; 

} 

if(DEBUG_ALIGN) fprintf(stdout, "ALIGN: largest piece match is %d monomers from 
%s\n" ,nlargestpiece,library [ilib] .name); 

if(DEBUG_ALIGN) fprintf(stdout/' ALIGN: Iargestpiece_ecount=%d, 
largestpiece_lcount==%d\n", 

largestpiece_ecount,largestpieceJcount); 

if(nlargestpiece >= smallest_acceptable_piece) { 

if(DEBUG_ALIGN) fprintf(stdout, tt ALIGN: incorporated^"); 
lcount = largestpiecejcount; 
ecount = largestpiecejxount; 
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while(ecount < (nlargestpiece + largestpiece_ecount)) { 

epothilone->alignedsequence[ecount] = 
libraryfilib] .monomersequenceflcount] ; 

strcpy(^othilone->alignedPKSname[ecoxint],library[ilib].name); 
strq>y(q)othilone->context[ecoimt],libiaiy 
epothilone->marked[ecount] =TRUE; 

if(ecount < (nlargestpiece + largestpiece_ecount - 1)) epothilone- 
>boundarytoright[ecount] = FALSE; 

lcount++; 
ecount++; 

> 

return (nlargestpiece); 

} /^maximal adj acentalignment*/ 

/* : 

PURPOSE: 

INPUT: 

OUTPUT: 

returns the size of the largest maximal adjacent set of monomers inserted. 
PROCEDURE: 

*/ 

int maximal_adjacent_alignment_and_dump( 
LIB *epothilone, 
int nwildcard, 

char wildcards[MAXWILD] [MAXLEN], 
LIB *library, 
int ilib, 

int smallest_acceptable_piece) 
{ 
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int ii=0,jj=0, kk-0; 

int ecount=0,lcount=0; 

int elen=0; 

int epothilonelen^O; 

int nlargestpiece=0,tnlargestpiece=0; 

int hold_this Jcount=0, hold_this_ecount=0; 

int wildcardmatch=FALSE; 

char *wptr; 

char *largestpiece_sptr,*largestpiece_eptr; 

char *hold_this ^lace^eptr, *hold_this_place_sptr; 

int Iargestpiece_lcount=0,largestpiece _ecount=0; 



char *sptr, *eptr, *lptr,*bufptr; 

ifi(DEBUG_WILDCARD) { 
iffawildcard > 0) { 

fprintf(stdout,"maximal_: wildcards[0]=%s\n M ,wildcards[0]); 

} 

fprintf^stdout/'maximal adjacent ^ 
smallest_accq}table__piece=%d\n T, ,smallest_acceptable_piece); 



sptr = library[ilib].monomersequence; 

eptr = epottalone->monomersequence; 

elen = strlen(epothilone->monomersequence); 

ecount=0; 

lcount=0; 

nlargestpiece=0; 

tnlargestpiece=0; 

hold_thisjplace_eptr = eptr; 

hold_this_ecount = ecount; 

while (*eptr !='\0'){ 

sptr = library[ilib].monomersequencc; 
Icount = 0; 
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/* NEW */ 
. hold_this_place_sptr = sptr; 
holdthisjcount = lcount; 
wildcardmatch = FALSE; 

while(*sptr != '\0 f ) { 

wildcardmatch = FALSE; 
if(epothilone->marked[ecount] — FALSE) { 

/* code for wildcards added MAS 05-16-00 */ 

wptr= 

if(*eptr = W) { wptr = wildcards[0]; } 
else if(*eptr === 'X 1 ) { wptr - wildcardsfl]; } 
else if(*eptr = Y 1 ) { wptr = wildcards[2]; } 
else if(*eptr = 'Z*) { wptr = wildcards[3]; } 

while(*wptr != W) { 

if(*wptr = *sptr){ 

wildcardmatch = TRUE; 
break; 

} 

wptr-H-; 

} 

if((wildcardmatch = TRUE) || (*eptr = *sptr)) { 
tnlargestpiece-H-; 

if(DEBUG_ALIGN) fprintf(stdout, "FOUND a match: 
len=%d, epo(%d, %c), lib[%d].name=%s (%d, %c)\n", 

talargestpiece, ecount, *eptr, 

ilib,library[i]ib].name,1couiit 3 *sptr); 

if(tnlargestpiece > nlargestpiece) { 

nlargestpiece = tnlargestpiece; 
largestpiecesptr = hold_this j>lace_sptr; 
largestpiece_Jcount = holdthisjcount; 
largestpiece_eptr = hold _thisj)Iace_eptr, 
largestpiece_ecount = hold this ecount; 
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if(DEBUG_ALIGN) fprintf(stdout, "FOUND a 
largest piece: len=%d, epo(%d, %c), lib(%d, %c)\n'\ 

nlargestpiece, largestpiece_ecount, 
*largestpiece_eptr. largestpiecelcount, *largestpiece_sptr); 

} 

sptr++; 
Icount-H-; 
eptrf+; 
ecount++; 

} else { 

tnlargestpiece = 0; 

sptrH-; 

lcount-H-; 

hold_this_p!ace_sptr — sptr; 
hold_this_lcount = Icount; 
/*NEW*/ 

eptr = hold_this_place_eptr; 
ecount = holdthisecoirat; 

} 

} else { 

tnlargestpiece = 0; 
break; 

> 

tnlargestpiece = 0; 
eptr = holdJhis_place_eptr + 1 ; 
ecount = hold_this_ecount + 1; 
hold_this_place_eptr = eptr; 
hold_this_ecount = ecount; 
if(DEBUG_ALIGN) { 

^rintf(stdout/ l iiicrementinghold - thisjplaceeptr=%c, 
hold^his_ecount^odW^*holdJhisjl^ 
} 



} 



/ 
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if(DEBUG_ALIGN) fprintf(stdout, n ALIGN: largest piece match is %d monomers from 
%s\n M ,nlargestpiece,library[ilib] .name); 

if(DEBUG_ALIGN) fprintffctdouV' ALIGN: largestpiece_ecount-%d, 
Iargestpiece_lcount=%d\n H , 

largestpiece_ecount,largestpiece_lcount); 

if(nlargestpiece >= smallest_acceptable_piece) { 

if(DEBUG_ALIGN) fprintf(stdouCALIGN: incoiporatedW); 

Icount = largestpiecelcount; 
ecount = largestpiece_ecount; 

fprintf(stdout, " ALIGN_TARGET: 
for(ii=0; ii<largestpiece_ecount; ii++) { 
Q>rintf(stdout, " "); 

} 

while(ecount < (nlargestpiece + largestpiece_ecount)) { 



epothilone->alignedsequence[ecount] = 
library[ilib].monomersequence[lcount]; 

strcpy(epothiIone->alignedPKSname[econnt],library[ilib].name); 
strcpy(epothilone->context[ecount],h"brary[ilib].context[lcount]); 
epothilone->marked[ ecount] = TRUE; 

if(ecount < (nlargestpiece + largestpiece ecount - 1)) epothilone- 
>boundarytoriglit[ecount] = FALSE; 

fprintf][stdout/ToC , \Ubrary[ilib].monomersequence[lcount]); 



lcount-H-; 
ecount-H-; 

} 



for(ii=ecount; ii<elen; ii++) { 
ftrintf(stdout," "); 
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} 

Q)rintf(stdout,"%s\n",library[ilib].name); 

} 

return (nlargestpiece); 
}/*maximal_adjacent_alignment_and_dump*/ 

/* 

PURPOSE: 
INPUT: 
OUTPUT: 
PROCEDURE: 

*/ 

int output_fresh_alignment( 

LIB *epotemp, 

int nboundary_cutoff) 

{ 

int acountK^ecount^O; 

int longest segmentlen=0,current segmentlen=0; 

char *aptr,*eptr; 

char boundary [MAXNAMELEN]; 



eptr = epotemp->monomerseqiience; 
ecount = 0; 

epotemp->nboundary = 0; 
longestsegmentlen = 0; 
currentsegmentlen = 0; 

strcpy(boundary,epotemp->alignedPKSname[ecovmt]); 
while(*eptr != W) { 

if(q)otemp->boundarytoright[ecount] = TRUE) { 
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q)Otemp->nboundar>HH-; 
if(current_segmentlen > longestsegmentlen) { 

longest_segmentlen = current_segmentlen; 

} 

currentsegmentlen = 0; 

} 

current_segmentlen++; 

ecount-H-; 

eptrH-; 

} 

if(current segmentlen > longest segmentlen) { 

longest_segmentlen = current_segmentlen; 

} 

if(epotemp->nboundary > nboundary_cutoff) return 1 ; 

eptr = epotemp->monomersequence; 
ecount = 0; 

fprintf(stdouV'HIT "); 
while(*eptr l=\0*) { 

if(epotemp->aJignedPKSname[ecount][0] — 1 W) { 

if(epotemp->bonndarytoright[ecount] — TRUE) { 

^rintf^stdout/'rocrTARGKyos)! ",*eptr,epotemp- 

>context[ecount]); 

} else { 

^rintflfstdout/^/ocrTARGC^os) ",*eptr,epotemp- 

>context[ecount]); 

} 

} else { 

if(epotemp->boundarytoright[ecount] = TRUE) { 

§)rintf(stdoiit; l %c:%4s(%s)[ ",*eptr,epotemp- 
>alignedPKSname[ecount],epotemp->context[ecount]); 
} e ^ se { 

^)rintf(stdout, ,, %c:%4s(%s) ",*eptr,epotemp- 
>alignedPKSname[ecount],q5otemp->context[ecoimt]); 
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} 

} 

ecount++; 
eptrH-; 

} 

jfcrmtf(stdout,"%d %d' , ,epotemp->nboimdary,longest__segmentlen); 
^rintfCstdout,""); 

eptr = epotemp->monomersequence; 
acount = 0; 
ecount = 0; 
while(*eptr != W) { 

if(epotemp->boundarytoright[ecount] — TRUE) { 

fprintf(stdout > "%c|",epotemp->context[acount][l]); 

} else { 

fprintf(stdout, M %c",epotemp->context[acount][l]); 

} 

eptrf+; 

acount-H-; 

ecount-H-; 

} 

fprintf(stdout/'\n"); 
return 1; 

} /*output_fresh_alignment*/ 

) 

/* 

*/ 

int get_library( 

char *libraryfile > 

LIB *library, 

int inter_modular_db_flag_off) 
{ 



r 
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int ii=Ojj=OJck=0,lcount=0; 

int nlib=0; 

char *bu^tr>uf[MAXBUF|; 

char tmonomersequence[MAXNAMELEN]; 

char *lptr,*tptr, 

FILE *libp; 

for(kk=0; kk < MAXJLIB_ENTRIES; kk++) { 
Ubrary[kk].recursion_tagged = FALSE; 
for(ii=0; ii<MAXNAMELEN; ii++) { 

libraxytkkJ.nametiiJ^W; 

library[kk].monomersequence[ii] = \0'; 

library[kk].alignedsequence[ii] = \&; 

for(ij=0; jj<MAXLEN; jj++) { 

library[kk].alignedPKSnameDj][ii] = '\0'; 
library[kk].markedOj] = FALSE; 
library[kk].context[0][O] = \0 f ; 
library[kk]xontext[jj][l] = W; 
hbrary[kk].context|jj]t2] = W; 
Ubrary[kk].context[jj][3] = W; 

} 

} 

library[kk].nboundary = 0; 

} 

/* read in the library from PKS.lib */ 
if(NULl^(libp=fopen(libraryfi { 

fprintf(stdout,"TRY AGAIN; couldn't open library file: %s\n",libraryfile); 
• nlib=0; 
exitO; 

} 

nlib=0; 

while(nlib < MAXJJB JENTRIES) { 
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iffNULL^getsCbuf^sizeof^uQ^ibp)) break; 
bufptr = buf; 

if(*bufptr = '#*) continue; 
if(*bufptr != '\n') { 

/* 

sscanf(bulptr, "%s %s 
%s",library [nlib] ,name,tmonomersequence ? library[nIib] .annotatedsequence); 

*/ 

Iptr = Iibrarytnlibj.name; 

while((*bufptr != ' ') && (*bufptr != W) && (*bufptr != V)){ 
*lptr++ = *bufptr++; 

} 

*lptr = W; 

if((*buQ>tr != \0') && (*bu$tr != V*)) bufptr++; 



lptr = library[nlib].monomersequence; 

while((*bufptr !=' ') && (*bufptr != W) && (*bufptr != V)) { 
/* 

This code specifically deletes inter-modular double bonds when the -d 

option is set. 

*/ 

, if(inter_modular_db_flag_off = TRUE){ 
if(*bufptr != '=•) { 

*lptr++ = *bufptr-H-; 

} else { 

bufptr++; 

} 

} else { 

*lptr++ = *bufptr++; 

) 

} 

*lptr = W; 

if((*burptr != \0') && (*bufi*r != V)) bufptr++; 
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lptr = library[nlib].annotatedsequence; 

while((*bufptr != ' •) && (*bufptr != \<W) && (*bufptr != V)) { 
*lptr++ = *bufptr++; 

} 

*lptr = W; 

if((*bufptr != \0') && (*bufptr != V)) bufptr++; 

Q)rintf(stdout 5 "LffiRARY(%(i) %s: 
%s\n"^ilib,library[nlib].name,library[nlib].monomersequence); 

fprintf(stdout,"LIBRARY(%d) %s: 
%s\n'\nlib,library[nlib].name,library[nlib].annotatedsequence); 



nlib++; 

} 

} 

fclose(]ibp); 

for(kk=0; kk< nlib; kk++) { 
lcount = 0; 

lptr = library [kk].monomersequence; 

while(*lptr!='\0'){ 

if(lcount = 0) { 

library[kk].context[lcount][0] = 
librarytkk].context[lcount][l] = *lptr; 
library[kk].context[lcount]f2] = *(3ptr + 1); 

} else { 

if(lcount = (strlen(library[kk].monomersequence) - !)){ 
Iibrary[kk]xontext[lcount][0] = *(lptr- 1); 
library[kk]xontext[lcount][l] = *lptr; 
library[kk].context[lcoxint][2] - »-'; 

} else { 

Ubrary[kk].context[lcount][0] = *0ptr - 1); 
library[kk]xontext[lcoiint][l] = *lptr; 
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library[kk].context[lcount][2] = *(lptr + 1); 

} 

} 

library[kk] .context [lcount] [3] = 'NO 1 ; , 

lptr++; 

lcount++; 

} 

fyrintf(stdout,"LIBRARY(%d) %s: 
%s\n",kk,library[kk] .name,library[kk] .monomersequence); 

^)rintf(stdout,"LIBRARY(%d) %s: 
%s\n",kk,library[kk].name,library[kk].annotatedsequence); 

for(jj=0; jj< strlen(library[kk].monomersequence); jj++) { 
Q>rintf(stdout,"(%s)",library[kk].context[jj]); 

} 

fprintf(stdout,"\n"); 

} 

return nlib; 
}/*get_library*/ 

/* 

*/ 

int dumpjSTARTER_align( 
LIB epothilone, 
int nwildcard, 

char wildcards[MAXWILD] [MAXLEN] 

) 

{ 

int elen=0; 

int ecount=OJtold_ecount=0; 
int wildcardmatch=FALSE; 
char *sptr,*eptr,*wptr; 

elen = strlen(epothilone.mononiersequence); 
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eptr = epothilone.monomersequence; 
fprintf(stdout, M ALIGN_TARGET: "); 
while(*eptr { 

§ointf(stdout > ,, %cVq)tr); 

epfaHH-; 

} 

^mtfCstdout,"^); 



ecount=0; 

eptr = epothilone.monomersequence; 
sptr = epothilone.alignedsequence; 
fprintf(stdout, "ALIGN_TARGET: "); 
while(*eptr != \V) { 

wildcardmatch = FALSE; 

wptr = ,m ; 

if(*eptr = 'X 1 ) { wptr = wildcards[0]; } 
else if(*eptr = V) { wptr = wildcards[l]; } 
else if(*eptr = 'Z') { wptr = wildcards[2]; } 



while(*wptr != T \0 f ) { 

if(*wptr=*sptr){ 

wildcardmatch = TRUE; 
break; 

wptrH-; 

} 

if(wildcardmatch = TRUE) { 
§>rintf(stdout, ,, : w ); 

} else { 

if(*eptr = *sptr) { 

fprintf{stdout, f T'); 

} else { 

fprintf(stdout, " "); 

} 

} 
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eptrf+; 
sptrH-; 

} 

fprintf(stdouV r Vn' r ); 

^rintfCstdout, "ALIGN_TARGET: "); 
eptr = epothiloncmonomersequence; 
ecount=0; 

while(*eptr != '\0') { 

if(epothilone.alignedsequence[ecount] = \0*) { 
fprintffctdout," "); 

} else { 

^rintf(stdout,"%c",epothilone.alignedsequence[ecount]); 
hold ecount = ecount; 

} 

eptr-H-; 
ecount-H-; 

} 

fprintf(stdout," %s\n H > epothilone.alignedPKSname[hoId_ecount]); 
fprint^stdout/STARTER^ALIGNiVa"); 

}/*dump_STARTER_align*/ 

/* 

*/ 

int reset_epotemp( 
LIB *epotemp 5 
LIB epothilone) 

{ 

int jj=0, elen=0; 

elen = sMen(epothilone.monomersequence); 
strcpyCepotemp-^ame^pothilone.name); 

strcpy(epotemp->monomersequence > q?othilone.nionoinersequence); 
strcpy(epotemp->aligaedsequence 3 epothilone.alignedsequence); 
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epotemp->nboundary = 0; 
for(ij=0; jj< elen; jj++) { 

strcpy(epotemp->alignedPKSname[jj],epothilone.aHgnedPKSname|jj]); 

epotemp->marked[jj] = epothilone.marked[jj]; 

strcpy(q>otemp->context[jj] 3 epothilone.context[jj]); 

epotemp->boundarytoright[jj] = epotWlone.boundarytoright[jj]; 

} 

} /*reset_epotemp*/ 

Thus, the present invention provides a useful means to generate new PKS genes and 
corresponding enzymes to produce polyketides. The invention having now been described by 
way of written description and examples, those of skill in the art will recognize that the 
invention can be practiced in a variety of embodiments and that the foregoing description and 
examples are for purposes of illustration and not limitation of the following claims. 

i 
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WHAT IS CLAIMED IS: 

1 . A method for representing the structure of a polyketide produced by a modular 
polyketide synthase, said method comprising the steps of: 

(a) defining a set of monomer units of which said polyketide is 
composed, 

(b) assigning an alphanumeric symbol or symbols to each different 
monomer unit in said set, 

(c) identifying one or more monomers in said set that is present in said 
polyketide, and 

(d) composing a string of said symbols ordered in a manner reflecting 
the order in which said monomers occurs in said polyketide, wherein said string 
of symbols represents the structure of said polyketide. 

2. The method of claim 1 , wherein said monomer set comprises two-carbon unit monomers, 
wherein a first carbon of said unit is substituted with hydrogen or methyl, and a second carbon 
of said unit is substituted with oxygen, hydroxy, or hydrogen, and said two carbon unit 
comprises either a single or a double bond between said first and second carbons. 

3. The method of claim 2, wherein said monomer set additionally comprises one or more 
members selected from the group consisting of two carbon unit monomers in which said first 
carbon is substituted with hydroxy, methoxy, or ethyl; a moiety corresponding to an amino acid 
or amino acid derivative incoiporated into a PKS by a non-ribosomal peptide synthase; a moiety 
corresponding to a structure incorporated into a polyketide by an AMP ligase or a CoA ligase; 
and a moiety corresponding to a structure corresponding to a structure in a polyketide after 
modification by a polyketide modification enzyme. 
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4. The method of claim 2 wherein the set of monomer unit and corresponding symbol 
comprises: 

OH OH OH OH 




or 



2 H OH o o 

=E ; J^ =F ' A Y /= G ; A^= H ; 



O 

A/ =| ; ^f=J; / Y= K ; -^ = l ; 




= M; and ^^ = N 



5. The method of claim 4 wherein the set of monomer unit further comprises a 
miscellaneous monomer that is assigned the symbol Q. 

6. The method of claim 4 wherein the set of monomer unit and corresponding symbol 
further comprises 

i 

2 H OH OH 



R 



OH 



R R 



= J,; = K ' ;and ""V = M " 

R ft U 
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7. A database of polyketides, in which each said member is represented by a string of 
alpha-numeric symbols, wherein said symbols represent structural subimits of said polyketide, 
and said string represents the order in which such subunits occur in said polyketide. 

8. The database of claim 7 that includes at least 1 00 different polyketides. 

9. The database of claim 7 wherein each said member is represented by a CHUCKLES 
string. 

1 0. The database of claim 7 wherein each said member is represented by an annotated 
CHUCKLES string. 

1 1 . The database of claim 7 wherein the symbol and its corresponding structural subunit are 
selected from the group consisting of 

OH OH OH . OH 




O 





and Q for a miscellaneous monomer. 
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12. The database of claim 7 wherein the symbol and its corresponding structural subunit are 
selected from the group consisting of 




oh oh 




and Q for a miscellaneous monomer. 

13. A database of polyketides, in which each said member is represented by a linearized 
representation of said polyketide. 

1 4. A method of designing a PKS gene capable of producing a desired polyketide, which 
method comprises: 

(a) defining a string of alphanumeric symbols representing the structure of said 
polyketide, 

(b) comparing said string to a database of strings of alphanumeric symbols 
representing polyketides produced by PKS genes, 
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(c) identifying common elements in said string representing the structure of said 
polyketide with elements in said strings in said database, and 

(d) generating one or more new strings from elements identified in step (b) that 
match said string representing the structure of said polyketide, wherein said new string defines a 
PKS gene capable of producing said polyketide. 

1 5. The method of claim 14, wherein all possible PKS genes encoding a desired polyketide 
from said database are generated and displayed. 

1 6. The method of claim 14, wherein said new strings generated in step (d) are rated and 
displayed in an order based on one or more parameters. 

1 7. The method of claim 16, wherein said parameters are selected from the group consisting 
of number of non-native module interfaces and number of non-native protein interfaces. 
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r 1= r 2 =H Erythromycin D 



Figure 1 
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Figure 2 
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CHUCKLES: ADGJDD 
SMILES: 

Cl(=0)-[C@H](C)[C@@H](OH)-IC@@H](C)[C@@H](0^ 

[C@@H](C)C-[C@@H](C)C(=0)-[C@H](C)[C@@H](C)- 

[C@@H](C)[C@@H](CC)01 



Figure 3 
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Figure 4A 
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Figure 4B 
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Figure 4C 



Design novel hybrid 
PKS genes from library / 
for TARGET / 




[initialize library counter2 to 1 ] 
\ set hokter_counter2 to counte^n ? 

t 



6il in ALtGNMENT blanks 
with longest stretches from 
entry based on Bbrary counter; ! 




increment 




library counter2 




4 YES 


if counter! equals countt 





Increment Bbrary 
counter^ 



write ALIGNMENT 
to file. 




copy aligned STARTER 
to ALIGNMENT again. 





increment library 




counter2 




setcounter2to 
holder_counter2 






increment 
hoWer_cout 


ater2 



This section allows the program 
to keep tiffing the ALIGNMENT 
until It Is completed 



This section gives equal weight 

to faring the blanks generated 

by the Initial STARTER alignment followed 

by any of the rest of me library 

components (traversed by counted). 



WO 01/92991 PCT/US01/1 7352 

7/9 



afign TARGET with 
STARTER and copy 
to ALIGNMENT 



allocate space for ALIGNMENT 
of same length as TARGET 



initialize TARGET ^pointer to start of TARGET 



initialize ALIGNLpointer to start of ALIGNMENT 



initialize a!igned_segmentjength to 0. 



initialize STARTER_pointer to start of STARTER 



"if TARGETjpointer a^>!£§ — ^SR**™ to mair i\ 
I of TARGET V^program 




increment align ad_segmenUenglh 



increment STARTER__pointer 



lfaB^n^_segmentJehgtii>= 1 
aligned_segmentjenb£h > 
ned^seg/nenUengJvfor ALIGNjiointer 



copy segment to ALIGNMENT 






copy alignecLse 
to ALIGNMENT 


gmenUength 






reset aligned_s< 


jgrnenLtength to 0 



increment TARGET_pointer 



Increment ALIGN_pointer 



Figure 4D 



8/9 



PCT/US01/17352 



Figure 4E 
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Figure 5 
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